LOCAL-IN-TIME EXISTENCE AND UNIQUENESS OF SOLUTIONS TO 
THE PRANDTL EQUATIONS BY ENERGY METHODS 



NADER MASMOUDI AND TAK KWONG WONG 

Abstract. We prove local existence and uniqueness for the two-dimensional Prandtl 
I ' system in weighted Sobolev spaces under the Oleinik's monotonicity assumption. In 

particular we do not use the Crocco transform. Our proof is based on a new nonlinear 
energy estimate for the Prandtl system. This new energy estimate is based on a cancellation 
property which is valid under the monotonicity assumption. To construct the solution, we 
use a regularization of the system that preserves this nonlinear structure. This new nonlinear 
structure may give some insight on the convergence properties from Navier-Stokes system 
to Euler system when the viscosity goes to zero. 
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1. Introduction 



The zero- viscosity limit of the incompressible Navier-Stokes system in a bounded domain, 
with Dirichlet boundary conditions, is one of the most challenging open problems in Fluid 
Mechanics. This is due to the formation of a boundary layer which appears because we can 
not impose the same Dirichlet boundary condition for the Euler equation. This boundary 
layer satisfies formally the Prandtl system. Indeed, in 1904, Prandtl [25] suggested that there 
exists a thin layer called boundary layer, where the solution u undergoes a sharp transition 
from a solution to the Euler system to the no-slip boundary condition u = on dfl of the 
Navier-Stokes system. In other words, Prandtl proved formally that the solution u of the 
Navier-Stokes system can be written as u = U + ubl where U solves the Euler system with 
U ■ n = on the boundary and ubl is small except near the boundary. In rescaled variables 
U + Ubl solves the Prandtl system. When studying this problem, there are at least 3 main 



(a) The local well-posedness of the Prandtl system; 

(b) Proving the convergence of solutions of the Navier-Stokes system towards a solution 
of the Euler system; 

(c) The justification of the boundary layer expansion. 

In full generality, these questions are still open except in the analytic case where (a)-(c) 
can be proved [23 ESI US EH] • 

Concerning (a), the main existence result is due to Oleinik who proved the local existence 
for the Prandtl system j22j |23] under a monotonicity assumption and using the Crocco 
transform (see also [3]). These solutions can be extended as global weak solutions if the 
pressure gradient is favorable (d x p < 0) [33j [31] • However, E and Engquist [6] proved a 
blow up result for the Prandtl system for some special type of initial data. More recently, 
Gerad-Varet and Dormy [7] proved ill-posedness for the linearized Prandtl equation around 
a nonmonotonic shear flow (see also [10J [8] ) . 

Concerning (b), the main result is a convergence criterion due to Kato [13] that basically 
says that convergence is equivalent to the fact that there is no dissipation in a very thin 
layer (of size v). This criterion was extended in different directions (see [291 EH E3 ) • Also, 



questions: 
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in [IT], it is proved that the convergence holds if the horizontal viscosity goes to zero slower 
than the vertical one. It is worth noting that the Prandtl system is the same in this case. 

Concerning (c), there is a negative result by Grenier [9] who proves that the expansion 
does not hold in W 1,co . Of course this does not prevent (b) from holding. 

There are many review papers about the inviscid limit of the Navier-Stokes in a bounded 
domain and the Prandtl system from different aspects (see (HEEDS]). Let us also mention 
that when considered in the whole space [281 IT2l 119] or with other boundary conditions such 
as Navier boundary condition j32j [TTJ [21 |2Q] or incoming flow [30], the convergence problem 
becomes simpler since there is no boundary layer or the boundary layer is stable. 

The prime objective of this paper is to prove the local existence and uniqueness for the two- 
dimensional Prandtl system under the Oleinik's monotonicity assumption in certain weighted 
energy spaces without using the Crocco transform. Precise statement will be provided in 
section |2j In addition to giving a very simple understanding of the monotonicity assumptions, 
our result may give us a better understanding about the questions (b) and (c) since it is 
given in physical space. Nevertheless, we are still not able to use our new nonlinear energy 
to study the convergence problem (b) or (c). In spirit this paper is similar to our previous 
paper about the Hydrostatic Euler equations [21] where we gave a proof of existence and 
uniqueness in physical space under some convexity assumption of the profile. The previous 
known proof of Brenier [3] uses Lagrangian coordinates and requires more assumptions on 
the initial data. 

Let us end this introduction by outlining the structure of this paper. In section |2] we will 
state our main result, that is, theorem 12.21 Explanations of our approach and approximate 
scheme will be provided in sections [3] and HI Assuming the solvability of approximate 
systems, we will derive our new weighted a priori estimates in section |5j Using these weighted 
estimates, we will complete the proof of our main theorem 12.21 in section [6j In section [TJ we 
will solve the approximate systems. For the sake of self-containedness, we will also provide 
several elementary proofs and computations in appendices \M- El Finally, let us mention that 
the new preprint p] also considers the existence for the Prandtl system in physical space. 
The methods of proof are very different. 



2. Main Result 



In this section we will first introduce the Prandtl equations, and then describe our solution 
spaces as well as our main result. Main difficulties and brief explanation of our approach 
will be given in section [31 
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Throughout this paper, we are concerned with the two-dimensional Prandtl equations in 
a periodic domain T x R + := {(x,y); x G R/Z, < y < +00}: 

' d t u + ud x u + vdyu = d 2 y u - d x p in [0, T] x T x R+ 

d x u + d y v = in [0, T] x T x R + 

(2.1) \ u\ t =o = uo on T x R + 

u\ y =o = v\ y= o = on [0, T] x T 

lim u(t, x,y) = U for all (t, x) G [0, T] x T, 

j/— >-+oo 

where the velocity field (u,v) := (u(t,x,y),v(t,x,y)) is an unknown, the initial data uq := 
u (x,y) and the outer flow U := U(t,x) are given and satisfy the compatibility conditions: 



(2.2) 



u L=o = and lim uq = U. 



Furthermore, the given scalar pressure p := p(t,x) and the outer flow U satisfy the well- 
known Bernoulli's law: 



(2.3) 



d t U + Ud x U 



-d x p. 



In this work, we will consider system (2.1) under the Oleinik's monotonicity assumption: 



(2.4) 



uj := d y u > 0. 



Under this hypothesis, one must further assume U > 0. 

Let us first introduce the function space in which the Prandtl equations (2.1) will be 
solved. Denoting the vorticity oj := d y u, we define the space H*'g for u by 



H s a 'J ■= J u : T x R+ -> R; \\cj\\ H ;i < +00, (1 + y) a u > 5 and ^ |(1 + y) a+aa D a u\ 2 < ^ 



\a\<2 



S 2 



where s > 4, 7 > 1, a > 7 + -, 8 G (0, 1), D a := d^d? 2 and the weighted H s norm || • 

2 y 

is defined by 
(2.5) 



l a, llff s .7(TxIR+) : ~ ~^ ?/) 7+a2 ^ )Q:u; llL 2 (Tx]R+)- 

a|<s 



Here, the main idea is adding an extra weight (1 + 2/) for each y- derivative. This corresponds 
to the weight — in the Hardy type inequality. Furthermore, we also denote H s,/y := {u : 

y 

T x R + -> R; IMInvr < +00}. 



Remark 2.1 (Requirement: cr > 7 + — ). If cr < 7 + — , then one may check that H*'g(T x 
is an empty set. Thus, we must have the hypothesis a > 7 + -. 
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Theorem 2.2 (Local H s a 'J Existence and Uniqueness to the Prandtl Equations (2.1)). Let 

s > 4 be an even integer, 7>1,<t>7 + - and 5 E (0, — ). For simplicity^, we suppose that 
the outer flow U satisfies 



(2.6) SUp |||Z7||| a +9,DO := SUp ||<9'f/|| H / s -2i+9,oo (T ) < +oo 



* 1=0 



Assume that uq — Ue H sn 1 and the initial vorticity ujq := d y uo E H s ^ 5 . In addition, when 
s = A, we further assume that 5 > is chosen small enough such that 

(2.7) Hlk" < cr 1 



where the norm || • ||# s ^ will be defined by (3.1) and C is a universal constant. Then there 



exist a time T := T(s,j,a,S, \\ujq\\h s ^,U) > and a unique classical solution (u,v) to the 
Prandtl equations plj] such that u-U E L°°([0, T]; H 3 ' 1 ^ 1 ) n C([0,T];H s - w) and the 



vorticity u := d y u E L°°([0,T]; H s a 'J) n C([0,T};H s - w), where H s - w is the space H s 
endowed with its weak topology. 

Remark 2.3 (U = constant). When the outer flow U is a constant, one may show that the 
life-span T stated in theorem 12.21 is independent of U. For the reasoning, see remark 16.31 

The proof of theorem 12.21 is based on our new weighted energy estimates, which relies on 
a nonlinear cancelation property that holds under the Oleinik's monotonicity assumption 



(2.4) An outline of our proof will be given in sections [3] and HJ and the detailed analysis will 
be provided in sections [5] - [71 

Before we proceed, let us comment on our notation. Throughout this paper, all constants 
C may be different in different lines. Subscript (s) of a constant illustrates the dependence 
of the constant, for example, C s is a constant depending on s only. 

3. Difficulties and Outline of Our Approach 

The aim of this section is to explain main difficulties of proving theorem 12.21 as well as our 
strategies for overcoming them. Let us begin by stating the main difficulties as follows. 



In order to solve the Prandtl equations (2.1) in certain H s spaces, we have to overcome 
the following three difficulties: 

(i) the vertical velocity v := —d~ x d x u creates a loss of x- derivative, so the standard 
energy estimates do not apply; 



1 The regularity hypothesis on the outer flow U is obviously not optimal in the viewpoint of our a priori 
weighted energy estimates. One may further loosen the regularity requirement on U by applying other 
approximate schemes. We leave this for the interested reader. 
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(ii) the unboundedness of the underlying physical domain TxE + allows certain quantities 
growth at y = +00 even if the solution is smooth or bounded in H s ; 

(iii) the lack of higher order boundary conditions at y = prevents us to apply the 
integration by parts in the variable y, but it is a standard and crucial step to deal 
with the operator dt — d%. 



Indeed, difficulty (i) is the major problem for the Prandtl equations (2.1) , and it explains 



why there are just a few existence results in the literature. The key ingredient of the current 



work is to develop a H s control by considering a special H s norm (see (3.1) below) which 
can avoid the regularity loss created by v. Difficulty (ii) is somewhat based on the fact 
that Poincare inequality does not hold for the unbounded domain T x M + . However, one 
may overcome this technical problem by replacing the Poincare type inequalities by Hardy 
type inequalities. This is our main reason for adding a weight (1 + y) for each ^/-derivative 



to our H s energies (2.5) and (3.1) Difficulty (iii) seems to be an obstacle, but it is not. 
A reconstruction argument for the higher order boundary conditions can fix this technical 
difficulty when s is even, see lemma 15.91 for more details. 

Now, let us explain our new weighted energy, which is the main novelty in this paper. 



Judging from nonlinear cancelations, the weighted norm (2.5) is not suitable for estimating 



solutions of the Prandtl equations (2.1) Thus, we introduce another weighted norm for the 
vorticity u, namely 



(3.1) 



where 



M 



ifJ'T(TxR+) 



;i+v) 7 *ii£. 



(TxR+) 



+ 



£ 

\a\<s 
ct\ <s— 1 



y y +a2 D a u\ 



L 2 (TxR+) 



diu-^-d s x {u-U) 



and 



u(t,x,y) := / u(t,x,y) dy 



provided that u := d y u > 0. The difference between norms (2.5) and (3.1) is that we replace 
the weighted L 2 norm of d^ui by that of g s , which is a better quantity because g s can avoid the 
loss of ^-derivative (i.e., difficulty (i) above), see subsubsection 15. L2l for further explanation. 
The first important observation is that as long as uj G H^.'], we can show that the new 



weighted norm (3.1) is almost equivalent to the weighed H s norm (2.5), that is, 



(3.2) 



MIh s < 7 ^ IM|h».t + \\ u - u\\ H s*->-i < IMIj#t + II^IU 2 



provided that u = d y u, u\ 



y=o 



and lim u 

y- 



U . The proof of (3.2) is elementary, and 



will be given in appendix lAl In spirit of (3.2), we will estimate ||u;||#j.7 instead of ||o;||^,7. 

The second important observation is that due to the nonlinear cancelation, the loss of x- 
derivative is avoided by the norm || • ||#|>7, so one can simply derive a priori energy estimates 
on uj by applying the standard energy methods. These estimates indeed can be extended to 
u e := d y u € , which is the regularized vorticity of the regularized Prandtl equations (4.1) below, 
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because the regularization (4.1) preserves the nonlinear structure of the original Prandtl 
equations (2.1) See subsection 15.11 for the detailed analysis. 

Once we have obtained the weighted energy estimates, it remains to derive weighted L°° 
controls on the lower order derivatives of u so that we can close our estimates in the function 
space H*'g. The derivations of these L°° estimates are standard: "viewing" the evolution 
equations of the lower order derivatives as "linear" parabolic equations with coefficients 
involving higher order terms that can be bounded by the weighted energies, we can obtain 
our desired estimates by the classical maximum principle since we have already controlled 
the weighted energies. These weighted L°° estimates are also extendable to the regularized 
vorticity u e , see subsection 15.21 for further details. 

In order to prove the existence, we will construct an approximate scheme which keeps the 
a priori estimates described above. Due to the nonlinear cancelation, our a priori estimates 
are complicated in certain sense, so the construction of the approximate scheme is tricky. 
An outline of this construction will be given in section HJ 

For the uniqueness, it is an immediate consequence of a L 2 comparison principle (see 
proposition 16. 4p . whose proof relies on a nonlinear cancelation that is similar to the one 
applied in the energy estimates. 



4. Approximate Scheme 



The main purpose of this section is outlining the approximate systems which we apply to 
prove the existence. Since our weighted H s a priori estimates are somewhat more nonlinear 
than usual, the approximate scheme is slightly more complicated. 

Our approximate scheme has three different levels and will be explained as follows. 

The first approximation of (2.1) is the regularized Prandtl equations: for any e > 0, 

' d t u e + u e d x u e + v e d y u e = e 2 d 2 x u e + d 2 u e - d x p e in [0, T] x T x R+ 



(4.1) 



d x u e + d v v 



;/ 1 
u%=o 



u 



y=0 



V 








lim u € (t, x,y) = U 

j/— !-+oo 

where p e and U satisfy a regularized Bernoulli's law: 

(4.2) d t U + Ud x U = e 2 d 2 x U - d x f. 



in [0, T\ x T x M+ 

on T x 1 + 

on [0,T] x T 

for all (t,x) G [0, T] x T, 



s 
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Or equivalently, the regularized vorticity co e : 
vorticity system: for any e > 0, 

{d t u e + u e d x u e + v e d y u e = e 2 d 2 x oo t + d 2 u e 
u% =0 = u) := dyU 
d y Uj e \y=Q = 8^ 

where the velocity field (u € ,v e ) is given by 



d y u e satisfies the following regularized 



in [0, T] x T x 
on T x R + 
on [0, T] x T 



+oo 



u e (t,x,y)dy and v e (t, x, y) := - / d x u e (t , x , y) dy 



(4.4) u^t,x,y) :=U 



The main idea of this approximation is adding the viscous terms e 2 dlu e and e 2 d 2 uj' L to avoid 
the loss of x-derivative. The advantage of this regularization is that our new weighted H s 
and L°° a priori estimates also hold for u € , and it is the main reason why we can derive the 

uniform (in e) estimates in section [51 The price that we pay is the appearance of extra terms 

d x oo e dlu e d x d v oo € , <9V , , „ , 
, , — and — - — during the estimation, but these terms can be controlled m 



a; 6 



uj'- 



oj'- 



the function space C([0, T]; H*'g). Before going to the next level, we should also emphasize 
that replacing the Bernoulli's law (2.3) by the regularized Bernoulli's law (4.2) is crucial 
here, otherwise the conditions u\ 



\y=o 



and lim u = U cannot be satisfied simultaneously. 

y— >+oo 



Although the approximate system (4.3) - (4.4) seems to be nice, its existence in the function 
space H^g is not obvious, so we will further approximate it by the next approximate system. 

The second level of approximation is the truncated and regularized vorticity system: for 
any e > and R > 1, 



d t u R + Xr{urO x uj r + v R d y u R } = e d x u R + d u R in [0, TjxTx 



(4.5) 



u R \ t =o = w := dyUQ on T x M + 

dyU R \ y=0 = d x p e on [0, T] x T 

where the velocity field (u R , v r ) is given by 



(4.6) u R (t,x,y):=U 



+oo 



u) R (t,x,y)dy and v R (t, x, y) := - / d x u R (t,x,y) dy 



Here, p e and U still satisfy the regularized Bernoulli's law (4.2) The cutoff function xr is 

y 

defined by Xr{v) := X\t>) wnere X ^ C™([0, +oo)) satisfies the following properties: 

R 



(4.7) 



0<X<1, xl[o,i] = l, suppx^ [0,2] and - 2 < X ' < 0. 



The main disadvantage of approximate system (4.5) - (4.6) is that the truncation on the 
convection term u R d x u R + v R d y ui R destroys the boundary condition w^l^o = as well as 
our weighted H s a priori estimate. However, it still keeps the weighted L°° controls. 
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To compensate for the lack of our new weighted H s estimates, one may apply the standard 



H s energy estimates because the system (4.5) - (4.6) does not have the problem of x-derivative 
loss. These estimates depend on e, but not on R. Thus, passing to the limit R — > +oo for 



the solution of (4.5) - (4.6) to that of (4.3) - (4.4) should be no doubt. The reason of doing 



this approximation is to prepare for our next approximate system. 

The third level of approximation is the linearized, truncated and regularized vorticity 
system: for any e > 0, R > 1 and n£N, 

r d t co n+1 + XR{u n d x u n + v n d y tu n } = e 2 d 2 x u n+1 + d 2 y co n+1 
(4.8) I Lu n+1 \ t=Q = u Q := d y u 

( d y u n+1 \ y=0 = d x p e 

where the velocity field (u n ,v n ) is given by 



in [0, T] x T x 
on T x R + 
on [0, T] x T 



(4.9) u n (t,x,y) :=U 



u n (t,x,y) dy and v n (t,x,y) 



d x u n (t,x,y) dy. 



In other words, (4.8) - (4.9) is a linearization of (4.5) - (4.6) 



The main advantage of the iterative scheme (4.8) - (4.9) is that its explicit solution formula 
can be obtained by the method of reflection. Using the explicit solution formula and the fact 
that XB,{ un dxU n + v n d y u n } has compact support, one may prove that there exists a uniform 
(in n) life-span T > for the approximate sequence {oj n } ne ^ provided that u G H*'^- This 
gives us a starting point so that we can solve the approximate systems and derive estimates. 

Solving the above approximate systems in a reverse order and deriving appropriate 
estimates, we can prove the existence to the Prandtl equations (2.1) Detailed analysis 



for solving the regularization (4.1) as well as other approximate systems (4.3) - (4.4), (4.5) 

(u e ,v 



(4.6) and (4.8) - (4.9) will be given in section [7J Assuming that 



,00 



solves (4.1 



(4.4) , we will derive uniform (in e) weighted estimates in section [5j Based on these uniform 



estimates, we will complete the proof of our main theorem 12.21 in section [6l 



5. Uniform Estimates on the Regularized Prandtl Equations 
In this section and the next we are going to complete the proof of our main theorem 12.21 



provided that we have a solution of the regularized Prandtl equations (4.1) In this section 



we will derive uniform estimates for the regularized Prandtl equations (4.1) by using the new 



weighted energy (3.1) introduced in section [3j These estimates are the main novelty of this 
paper. Then we will finish the proof of theorem 12.21 in section [61 After that, an outline for 



solving the regularized Prandtl equations (4.1) will be provided in section [7J 

Our starting point is that we can solve the velocity field (u e ,v e ) from the regularized 



Prandtl equations (4.1) More precisely, let us assume proposition 15.11 which will be shown 



in section [7J below for the moment. 
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Proposition 5.1 (Local Existence of the Regularized Prandtl Equations). Let s > 4 be 

11 

an even integer, 7 > l,er > 7 + -,6 G (0, -) and e G (0,1]. If ujq G H s a2S ' 7 , U and 



p e are given and satisfy the regularized Bernoulli's law (4.2) and the regularity assumption 



(2.6), then there exist a time T := T(s, 7, a, 5, e, ||u; ||#s+4, 7 , U) > and a solution uj e G 
C([0,T];#*"£ 4,7 ) n C\[0,T];H s+2 ' J ) to the regularized vorticity system pgWO) 



Furthermore, the velocity (u e ,v e ) defined by (4.4) satisfies the regularized Prandtl equations 



(4.1) as well. 



Remark 5.2 (Initial Data). The H^'^ functions can be approximated by H s a + ^ n functions 
in the norm || ■ ||#s, 7 , so by the standard density argument, the hypothesis ojq G H^s' 2 '' 7 can 
be reduced to be ojq G H s j1 & in our final result. 



According to proposition 15.11 the life-span T S)T)0 - l( s )eM)l [/ of u e depends on e, so our aim in 
this section is to remove the e-dependence by deriving uniform (in e) estimates on a/. In 
other words, we will prove the following 



Proposition 5.3 (Uniform Estimates on the Regularized Prandtl Equations). Let s > 4 

1 



be an even integer, 7 > 1, a > 7 + -, 5 G (0,1) and e G [0,1]. If uj e G C([Q,T};H^. 



n 



C 1 ([0, T]; H s a +2a ) and (u e ,v e ,u £ ) solves (4.1) - (4.4), then we have the following uniform (in 
e) estimates: 

(i) (Weighted H s Estimates) 



(5.1) 



< <\\uJo\\ H ^+ / F(T)dr} il-C mS [\\oj / F(T)dT) t 



8-2 
2 



1 

" s-2 



as long as the second braces on the right hand side of (5.1) is positive, where the 
positive constant C Sj1j(Tj s depends on s, 7, a, 5 only and F : [0,T] — > IR + is defined by 



(5.2) 



F := C., W {1 + ll^f/HU + C s ll^|| 2 H 3- 2i(T) . 



(ii) (Weighted L°° Estimates) Define I(t) := ^ |(1 + y) a+a2 D a u £ (t)\ 2 . For any s > A, 

H<2 



(5.3) 



|/(f)IU~(TxR + ) < max{||/(0)|| L co (TxR+) ,6C 2 fi(t) 2 }e c ^-^ 1+G W> t 



where the universal constant C is the same as the one in inequality (B.3) , Q and 
G : [0, T\ — > R + are defined by 

(5.4) Q(t) := sup ||w e ||_H- s ^ and G{t) : = sup ||a/||# s ,7 + sup ||9|[/||i2( T ). 

[o,t] " [o,t] " [o,t] 



(5.5) 



(5.6) 
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In addition, if s > 6, then we also have 

||/(t)IU«crxH+) < (||/(0)|| L oo (TxR+) + L7 Si7 {l + fi(t)}fi(t) 2 t)e c --^ 1+G W^. 
For s > A, we have the following lower bound estimate: 



n 



min(l + yYujHt) 

TxR+ 

> (1-L7 s , w {l + C7(t)}te 



C s ^^ s {l+G(t)}t 



) • f min (1 + y) a u - C S)7 n(t)t ) 



provided that min(l + y) a co — C sa Q(t)t > 0, where C s>J>CTi s is a positive constant 
depending on s, 7, a and 5 only. 



Remark 5.4 (Two L°° Estimates on I). In proposition 15.31 we stated two L°° controls on the 



quantity I(t), namely, estimates (5.3) and (5.5) Indeed, (5.5) is a better estimate within 



a short time, but it only holds for s > 6. Thanks to this better estimate, we can derive 



the uniform weighted L°° bound (6.2) without any additional assumption when s > 6. In 



contrast, we are required to impose an extra initial hypothesis (2.7) for the case s = 4 since 



we only have the weaker estimate (5.3) in this case. See proposition 16. II for the details. 



Remark 5.5 (A Priori Estimates on the Prandtl Equations). When e = 0, proposition 15.31 



provides a priori estimates for the Prandtl equations (2.1) Similar situation occurs in 
proposition 16.11 as well. 



The proof of proposition 15.31 will be given in the subsections 15.11 and 15.21 as follows. 



5.1. Weighted Energy Estimates. The objective of this subsection is to derive uniform 
(in e) weighted H s estimates on u e . These estimates, which are the main novelty of this paper, 
include: (i) weighted L 2 estimates on D a u € for \a\ < s and a\ < s — 1 in subsubsection 
15.1.11 and (ii) weighted L 2 estimate on g s in subsubsection 15.1.21 We will combine these two 



estimates in subsubsection l5~TT3l to obtain the uniform weighted energy estimates (5.1) This 
will complete the proof of part (i) of proposition 15.31 



5.1.1. Weighted L 2 Estimates on D a u € . Using the standard energy method, we will derive 
weighted L 2 estimates on D a u e for \a\ < s and ai < s — 1 in this subsubsection. It works 
because we are allowed to loss at least one x-regularity in these cases. 
More specifically, we will prove 
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Proposition 5.6 (L 2 Controls on (1 + y) J+a2 D a co e for \a\ < s and ct\ < s — 1). Under the 
hypotheses of proposition 13731 we have the following estimates: 

~ E \\(i+yr a2 D^\\h 

\a\<s 
«i<s- 1 

(5.7) H< s M< s 

v ' ai<s— 1 ai<s— 1 



.e 1 1 2 

*9 



+ C s>7><7iS {\\u e \\ H sn + \\d s x U\\ L ^}\\oj e \\jj 

+ C S , W {1 + ||a; e ||^7} s - 2 ||a; e ||^, 7 + C s ^ II^P' 
where the positive constants C s and C S)1)(T ^ are independent of e. 



ell 2 

«=0 



Remark 5.7 (Boundary Terms at y = +oo). In the proof of proposition 15.61 and that 
of proposition 15.101 below, we will ignore the boundary terms at y — +oo while we are 
integrating by parts in the variable y. Skipping these boundary terms is just for the 
presentation convenience, and ignoring these technicalities is harmless. Indeed, one may 
deal with these boundary terms by any one of the following two methods: 



(i) Since u> e G H s a + 5 in , by proposition IC.1[ we have nice pointwise decays (C.l) for oj € 
and its spatial derivatives. Therefore, when a is much larger than 7, one may easily 
check that those terms which we will omit actually vanish; 

(ii) As long as oj e (t) G Hp], the norm ||o; <E ||^.7 < +00 provides certain integrability of 
the underlying quantities. Thus, one may overcome the technical difficulty by first 

y 

multiplying by a nice cutoff function Xr{v) '■= x{~^) during the estimation, and then 

R 

passing to the limit R — > +00. The main advantage of this approach is that it only 
requires a > 7 + - and w e (t) in Hp], but not in H^'" 1 . As a demonstration, we will 
apply this argument in the proof of proposition 16.41 for the reader's convenience. 

In conclusion, the proofs of proposition 15.61 and 15. 101 are absolutely correct, even if we ignore 
the boundary terms at y = +00. 



Proof of proposition 15. 6\ Differentiating the vorticity equation (4.3) L with respect to x ot\ 
times and y CK2 times, we obtain the evolution equation for D a co e : 



(5.8) {d t +u e d x +v e d y -e 2 d 2 x -d 2 v }D a uj 



O<0<a 
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Multiplying |(5.8)| by (1 + y) 2l+2a2 D a u e , and then integrating over T x M + , we have 



lj t m+yy +a2 D a ^\\i 

e 2 [[(1 + y) Wa2 D a u e d 2 x D a u t + / / (1 + y) Wa2 D a u e d 2 D 



(5.9) 



1 + y)^ +2a2 D a u € {u e d x D a u e + v e d y D a u e } 



- [a] 1 7 \l + y) 2ry+2a2 D a u e {D^u e d x D a -^u e + D^v e d y D a -^u e }. 
Now that we can apply integration by parts and the standard Sobolev's type estimates on 



trilinear forms to control the right hand side of (5.9) as follows 



Claim 5.8. There exist constants C sn and C arft(T £ > such that for any \a\ < s and 
a i < s — 1, 



(5.10) 



(5.H) 



(5.12) 



(5.13) 



e 2 / / (1 + y) Wa2 D a u e d 2 x D a u e = -e 2 ||(l + y y +a2 d x D a u e \\ 2 



L 2 • 



JJ0- + y) 2l+2a2 D a co e d 2 y D a Lo e 

< --\\{l+yY <+a2 d y D a u t \\ 2 L 2- [ D a u e d y D a u e dx\ =0 + C sn \\u e \\ 2 H s,. 

JJ{1 + y) Wa2 D a u e {u e d x D a uj e + v e d y D a oj e } 
< C sa><7iS {\\u*\\ Hr + \\d:U\\ L z}\\u*\\ 2 Hr 

J2 jj (1 + y) 2l+2a2 D a u t {D^u t d x D a -^u t + D^v e d y D a -^u e } 



0</3<ct 



Assuming claim I5.8[ which will be shown later in this section, for the moment, we can 



apply inequalities (5.10) - (5.13) to the equality (5.9), and obtain 
(5.14) 

lj t m+yy +a2 D^\\h 

< -e 2 \\(l + yy +a2 d x D a ^\\ L 2-^\\(l + yy +a2 d y D a u*\\ 2 L2 



^D a u%D a ^ dx\ y=0 + C Sl7 IM|^,, + C s ^ 5 {\\u*\\ Hr + ||^|| Loo }||a; £ || 2 Hr . 
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When | a | < s — 1, we can apply the simple trace estimate 



(5.15) 



|/| dx 



< C 



y=o 



JT 



|/| dxdy + 



\d y f\ dxdy 



to control the boundary integral / D a tu e d y D a uj e dx 

1 



JT 

as follows: 



(5.16) 



D a u e d y D a u e dx\ y=0 



12" r ' *' y 

However, when |a| = s, a main difficulty arises: the order of d y D a u e \ y= o is too high so that 

we cannot control the boundary integral / D a u e d y D a u e dx\ y= o by the simple trace estimate 
./ 



(5.15) In order to make use of (5.15), we must reduce the order of the problematic term 
d y D a u e \ y= Q. When s is even, this can be done by a boundary reduction argument as follows. 

At this moment, let us state without proof the following boundary reduction lemma, which 
will be proven at the end of this subsubsection. 



Lemma 5.9 (Reduction of Boundary Data). Under the hypotheses of proposition HOI we 
have at the boundary y = 0, 

d y uj e \ y=0 = d x p e 



(5.17) 



d 3 y uj% =0 = {d t - t 2 d 2 x )d x f + LO e d x u% =0 . 



For any 2 < k < [—], there are some constants C k i p \ ( p. ... ^ 's, which do not depend on e or 
(u e ,v e ,u e ), such that 

k—l max{2,fc— 1} j 

(5.18) df +1 co% =0 = {dt - e 2 dl) k d x f + Y, e21 E C ^uP,- V ft^Vo 

1=0 j=2 p£ A l,i i=1 
j j j 



where A j k l := {p := (p 1 , p 2 , • • • , p J ) G N 2j ; 3 ^ p\ + ]T p\ = 2k + 41 + 1, ]T p\ < k + 21 - I, 

i=i i=i i=i 

j 

^2pl < 2k -21 - 2 and \p l \ < 2k - I - 1 for all i = 1,2,..., j}. 



i=i 



Now, we can apply lemma [5791 to control the boundary integral / D a uj e d y D a uj e dx\ y=0 for 
| a | = s with < a± < s — 1 in the following two cases: 



Case I: (c*2 is even) 
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When a 2 := 2k for some k G N, we can apply boundary reduction lemma 15.91 to 
y=0 , and obtain 



d y D a u 



D a u e d y D a u) £ dx\ y=0 



(5.19) 



J D a ^(d t -e 2 d 2 x ) h d^ +1 fdx\ y=0 



fe_l max{2,fc-Z} / j \ 

+ E e2 ' E E ^^.-v / ri DPV ^=°- 



1=0 



According to the definition of A 3 kl , one may check by using the indices restrictions that the 
largest possible order for d" 1 D p uj e is < s — 1 and at most one of d" 1 D p uj e can attain the 
order s — 1, namely, the orders of other terms are < s — 2. Therefore, we can apply the 



simple trace estimate (5.15) and proposition IB. 31 to the identity (5.19) to obtain 
(5.20) 



D a u e d y D a uj e dx 



y=0 



< 



^H(i + y y +a >d y D a ^\\ 2 L2 + c s J2 \\%W B -*<x) + 



I, ,e|| \s-2||. ,e\\2 



Case II: (a 2 is odd) 

When a?2 := 2k + 1 for some k G N, since a.\ + a 2 = s is assumed to be even, we know 
that ai > 1. Using integration by parts in x, we have 



(5.21) 



D a uj e d y D a uj e dx 



d r D a uW 1 d^ +1 u t dx 



y=o 



Now, the term d x D a u e \ y =o = ^ e \y=o has an odd number of y derivatives, and 

hence, we can apply the boundary reduction lemma I5\9l to reduce the order of the right hand 



side of (5.21) Similar to the Case I, we can further apply the simple trace estimates (5.15) 



and proposition IB. 31 to eventually obtain the following estimates: 
(5.22) 



D a u e d y D a u e dx 



y=o 



< ^ll(i + y)^ +1 ^ 1_1 ^ + V||i a + ^^11^11^^ 

1=0 



s-2||, ,e\\2 

g J il \\H„ 



Finally, combining estimates (5.14), (5.16), (5.20) and (5.22) and summing over a, we 

□ 



prove (5.7) 
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In order to complete the proof of proposition 15.61 it remains to show claim 15.81 and the 
boundary reduction lemma 15.91 Let us first prove the claim 15.81 as follows. 



Proof of claim \5.8\ . 



Proof of (5.10) 



The equality (5.10) follows immediately from an integration by parts in the variable x. 



Proof of [(531) 

Integrating by parts in y (cf. remark I5TT1) . we have 

J J {I + y) Wa2 D a u e d 2 D a u e 

= - \\(l + yy +a2 d y D a u e \\ 2 L2 - [ D a co e d y D a u e dx\ 

J 



y=0 



- 2( 7 + a 2 ) / / (1 + y)^ +2a '- 1 D a uj e d y D a 



to 



< -~\\(l + y y +a *d y D a u e \\l 2 - [ D a ^d y D a u e dx\ y=0 + C s 
4 Jj 



which is inequality (5.11 



Proof of (5.12) 



Integrating by parts (cf. remark ISTTI) . and using d x u e + d y v e = 0, we have 
(5.23) 

(1 + yf' <+2a2 D a uj e {u t d x D a uj e + v e d y D a co e } = (7 + a 2 ) [ [ (I + yf^^v^D 



which and inequality (B.9) imply inequality (5.12) 



Proof of (5.13) 



Using the facts that d y v e = —d x u e and d y u e = u e , one may check that all terms on the left 



hand side of (5.13) are one of the following three types: denoting ei := (1, 0) and e 2 := (0, 1), 
for 77 G N and k,6 G N 2 , 

Type I: 

j x := JJ{^ + y) 2l+2a2 D a u e d v x v e D K tu e 
where 1 < r\ < s — 1 and rjei + k = a + 



Type II: 



j 2 := JJ(l + y) 2j+2a2 D a u e d v x u e D K u e 
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Type III: 



Ja := 



1 + y) 2 ~ t+2a2 D a u € D e u € D K t 



where \9\ < s — 1 and 6 + k = a + e\ — e 2 . 



Thus, it suffices to control Ji, J 2 and J 3 by the right hand side of (5.13) as follows. 



Estimates for Type I: 

When 1 < 7] < s — 2, applying proposition IB.3[ we have, since k 2 = a 2 + 1, 



\-h\ < \\(i + y y +a2 D a ^\\ L 2 



l + y 



[l+y)^D^\\ L , 



< C s , w {||a;1^ + ||^|| L2 }||a;1^ 



When rj — s — 1, by triangle inequality and proposition IB.3| we have, since k 2 — ot 2 + 1, 



< \\(l + yV +a2 D a ^\\ L , 



d s -^ + yd s x U 



l + y 



{l+yy +K2 D K u e \\ LO 



L- 



+ n(i + ^ +a2 wi| L2 wml- \\{i + y y + ^D^\\ L , 



In conclusion, J\ can be controlled by the right hand side of (5.13) 



Estimates for Type II: 

When 1 < 77 < s — 1, applying proposition IB.3[ we have, since n 2 = «2? 



|J 2 | < + yV +a2 D a ^\\ L 2\\d2u e \\ L ^\\(l + y y +K w K ^\\ L 2 



< ^,,{11^11^ + ||^£7|| La }||a; 



e II 2 
I A3 



When 77 = s, by triangle inequality and proposition IB.3[ we have, since k = (0, a 2 ) and 

< a 2 < 1, 

|J 2 | < ||(l + y) 7 ^^a; e ||^||^(w-C0||^||(l + y) 7+ ^a; e |U- 
+ ||(l + ^ +Q2 ^l L2 ||^f/|| L o ||(l + 2/ )^ J D K w l L2 
< a i7i<J|5 {||o; e ||^,-r + ||^C/|| ioo }||a; e ||^. 



In conclusion, J 2 can also be controlled by the right hand side of (5.13) 



Estimates for Type III: 
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Applying proposition IB. 3\ we have, since 9 2 + k 2 = a 2 — 1, 
< CwAWu'WeT + \\d s x U\\ L2 }\\u% s , 7 . 



Combining all estimates for type I - III, we prove inequality (5.13) 



if |0| < s -2 
if 101 = s - 1 



□ 



Lastly, we will prove the boundary reduction lemma 15.91 as follows. 



Proof of lemma UrU First of all, let us mention that equality (5.17) L is exactly the same as 
the given boundary condition (4.3) y Furthermore, differentiating the vorticity (4.3) L with 
respect to y, and then evaluating at y = 0, we obtain equality (5.17) 2 by using (5.17) L and 



\y=o 



\y=0 



0. Thus, it remains to prove the formula (5.18^ 



In order to illustrate the idea, let us derive the formula (5.18) for the case k = 2 as follows. 
Differentiating the vorticity equation (4.3) L with respect to y thrice, and then evaluating 

■; e |y=o = 0, 



at y — 0, we obtain, by using (5.17) } and u e \ y= o 
(5.24) 

l»=o 



<9VL=n = (d t - e 2 d 2 x ) 2 d x p e + (d t - e 2 d 2 x )(^d x ^) + 3u e d x d 2 v u e + 2d y oj e d x d y u e - 2d x u e d 2 u e \ 



y=o- 



Since the last three terms on the right hand side are our desired forms, we only need to deal 
with the terms (dt — e 2 (9^)(o; e 9 x a; e )| l/= o. Using the evolution equations for co e and d x co e as 



well as u e 
(5.25) 



y=0 



\y=o 



= 0, one may check that 



(d t -e 2 d 2 x )(^d x ^ 



\y=0 



u e d x d 2 uj e + d x u e d z tu e - 2e 2 d x u e d 2 x u e 



y=o, 



where all terms on the right hand side of (5.25) are also our desired forms. Substituting 
(5.25) into (5.24), we justify the formula (5.18) for k = 2. 



Now, using the same algorithm, we are going to prove the formula (5.18) by induction on 



k. For the notational convenience, we denote 



k-l 



max{2,fc— /} 



k,l,p l ,p 2 ,- ,pi 



to'- 



=o}. 



1=0 



i=2 



1=1 



Under this notation, we will prove d 2k+l uj e \ y= Q — (d t — e 2 d x ) k d x p e G Ak- 

Assuming that the formula (5.18) holds for k = n, we will show that it also holds for 
k = n + 1 as follows. 
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In order to reduce the order of d 2n+3 ui e \ y= o, we first differentiate the vorticity (4.3) L with 
respect to y 2n + 1 times, and then evaluate the resulting equation at y = to obtain 

2rc+l 



(5.26) 



d^u% =0 = (d t - e 2 d 2 )^+V + £ ( 

m=l ^ 

-E 



2n + 1 



di™-V&.a 2ri - m+ V 



m=2 



2n + 1 
m 



<9 x .d™-V«9f- m+ V| 



y=o- 



By a routine checking, one may show that the last two terms of (5.26) belong to A n+ i, so it 
remains to deal with (d t — e 2 dl)d 2n+1 u e \ y=0 only. 



Thanks to the induction hypothesis, there exist constants C nj i tP i jP 2 t ... >p} -'s such that 

n— 1 max{2,n-£} j 

d^u;% =0 = (d t - e 2 dlYd x f + E^ E E C n,i,P^,-,^ 

1=2 p^K, i=1 



1=0 



so we have, up to a relabeling of the indices p J 's, 
(d t -e 2 dl)d? +1 uj% =0 

n—l max{2,n— /} j 

= (d t - e 2 dir^ +J2e 21 Yl E r >■■>;■ -^ & - ^iw 1 ^ n Dpl 

(5.27) i=o 3=2 peA i 



i=2 



-1 max{2,ra— 1} 



-E 

z=o 



.2Z+2 



E E ^rM,P 1 ,P 2 ,- >/>> 

i=2 pe<, 



d x D pl oo e d x D p2 uj t f[D pt 



i=3 



where C n ^^ p i^2 ... p j's and C„^ p i iP 2 r . ^'s are some new constants depending on C n ^ i pi i p2 j ... p i's. 
It is worth to note that the last term on the right hand side of (5.27) belongs to A n+ i, so it 
remains to check whether the second term on right hand side of (5.27) also belongs to *4 n+ i. 
Differentiating the vorticity equation (4.3) L with respect to x p\ times and y p\ times, and 



then evaluating at y — 0, we have, by using u 
(5.28) 

fa ,2x2\ no 1 . .£ 

\y=0 



\y=o 



v e \ y= o = and denoting e 2 '■= (0, 1), 



{dt - t 2 d 2 x )D pl uj e \ 

- (f^D^^d x D pl -?^ + ( pl )d x D^u%D p1 ^^ + d 2 D pl co% =0 . 



^P 1 
&>1 



^p 1 

/3 2 >2 



Using (5.28) , one may justify by a routine counting of indices that the second term on the 

□ 



right hand side of (5.27) belongs to A n+ i. This completes the proof of lemma 1531 
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5.1.2. Weighted L 2 Estimate on g e s . In this subsubsection we will derive the L 2 estimate 
on (1 + y^gl by the standard energy method. This can be done because the quantity 

U) avoids the loss of x-derivative by a nonlinear cancellation, 



91 



diu e 



— — d x (u e 



which is one of the key observations in this paper and will be explained as follows. 
Let us begin by writing down the evolution equations for u e and u e — U: 



(5.29) 



(d t + u e d x + v e d y )u e = e 2 d 2 u € + <9, 
{d t + u e d x + v e d y )(u e -U) = e 2 d 2 x (u e -U) + d 2 (u e - U) 



u 



U)d x U 



where we applied the regularized Bernoulli's law (4.2) in the derivation of (5.29) Since our 



aim is to control the H s norm of u e (or u e — U), let us differentiate (5.29) with respect to 
x s times. Then we have 



(5.30) 



{dt + u e d x + v e d y )d s x cj e + d s x v e d y u e = e 2 d s x +2 u e + d s x d 2 co e + 



(d t + u € d x + v e d y )d s x (u € -U) + d s x v e u e = e 2 d: 



•+2, 



where we applied the fact that uj 6 = d y (u e — U) and the symbol 
order terms which we want the reader to ignore at this moment. 



U) + d s x ^-U) + --- , 
represents the lower 



The main obstacle in (5.30) is the term d x v e = —d x d~ u e } which has s + 1 ^-derivatives 
so that standard energy estimates cannot be closed. However, since there are two equations 
in 



(5.30) we can eliminate the problematic term d x v e by subtracting them in an appropriate 



way. 



Subtracting 



(5.30) , from (5.30) L , we obtain 



(5.31) 



(d t + u e d x + v%-e 2 d 2 x -d 2 y )gl 

2e 2 K + V - CO - - C/)Ra e + 2^ - g\d s x U - £ f'W^" 

j= i v/ 



CO' 



s-l 



3=1 



s-l 



E .)d s - j v e {did y u* - a'<)^<} ■ WiV - U)d s x ^ +1 U. 



3=0 



where g% := d x u e - a e d x (u e - U) and a e : 



d v u e 



ijj'- 



. Here, the main reason that we can apply 



this nonlinear cancelation is the Oleinik's monotonicity assumption (i.e., u > 0), which is 



ensured in our solution class H^J. For the justification of (5.31), see appendix [Dl 
Now, we are going to derive the following weighted energy estimate on g e s : 
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Proposition 5.10 (L 2 Control on (1 + yyg e s ). Under the hypotheses of proposition HOI we 
have the following estimate: 

(5.32) 



2dt' 



< - ^||(i + yyd x9 !\\h - + yyd y gl\\h + cwarvwhm + c^WMi 



ei|2 



= (T) 



e || 2 



+ C s>1i<7i s{1 + \\oj e \\ H ^ + \\d s x U\\ L ~ m }{\\u e \\ Hr + ||^ +1 t/||L«cr)} IMIhj 
where the positive constants C,C j7 s and C arij<Tt $ are independent of e. 

The proof of proposition 15.101 is almost a straight forward application of energy methods 



except the estimation (5.37) below is slightly tricky. 

Proof of proposition UTTU . Multiplying the evolution equation (5.31) by (l + y) 2y g e s , and then 
integrating over T x R + , we have 

(5.33) 



2dt 

2 



JJ(1 + yf'gldlgl + //(! + y^g^gt — J J (1 + y)*><f.{u<d ai gl + v%gl} 



+ 2e 2 I l(l + y)^gl{d s x +1 (^-U) 



LO'- 



d s x (u e -U)}d x a £ 



+ 2 + y)^\gl\ 2 d y a^ - JJ{1 + y)^g{gldlU 



s-l 



£ Q //(! + y) 2l 9>^ - U)d^ +1 U. 



Indeed, all terms on the right hand side of (5.33) can be controlled by using integration 
by parts and the standard Sobolev's type estimate on multilinear forms. Precisely, we have 

Claim 5.11. There exist constants C,Cg,C 7t g and C Si7iCTj 5 > such that 
(5.34) 



/ 1 (i + y) 2 " 'g&st = -t 2 \\(i + yyd x9 :\\h, 



(5.35) 



-1-1 e II 2 

r P \\L 2 (J)l 



(5.36) 



JJ(l + y)^gl{u^ x gl + v%gl} 
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(5.37) 



(5.38) 



(5.39) 



(5.40) 



(5.41) 



2e 2 11(1 + 2/) 27 £|R + V - tO - ^d S M - U)}d x a 



d r u) e 



id 1 - 



< -6 2 ||(1 + yyd x gl\\l 2 + e 2 C s ^ 5 {\\^\\H^ + ll^ +1 C/||x 2 (T)}|k e ||^7, 



2 //(l + y)W 



<<^IMIk", 



|y (1 + v)^g\gtdlU 



<c s ,,, a 4d s x uh- m m\\ Hr + 



< C SjW {||o; e || H| „ + Rl7||L» m rlMlflr 



< 



(5.42) 



s-l 

E 

3=0 



CwflMlnT + ||^||L<-(T)}||w e ||^, 



{l + y)^glcfdi{u*-U)d s -i +x U 



< C s , w ||^ +1 ^|U» (T) {||a; e || Hr + ||^C/|U 2(T) }||a; e |U r . 
Assuming claim |5~TT| which will be proven at the end of this subsection, for the moment, 



we can apply (5.34) - (5.42) to (5.33) to obtain our desired inequality (5.32) because e G [0, 1] 

and \\%u\\v<n < WMl-vo < ll^l^m ^ ll^ +lf/ IU-m- 



□ 



To complete the proof of proposition I5.10[ we will show claim 15.111 as follows. 
Proof of Claim \5.11\ 



Proof of (5.34) 



The equality (5.34) follows directly from an integration by parts in the variable x. 



Proof of (5.35) 



Integrating by parts in the variable y (cf. remark [5T7|) . we have 
JJ(l + y) 2l 9ld 2 y gl 
= -\\(l+ y y dy gl\\l 2 - ^gtd y gtdx\ y=0 -2 1 j j (1 + y) 2ry ~ l g\d y gl 

(5.43) < — ^||C1 -H 2/)^^^!]^ — ^ S'S^/S'S ^1^=0 CVM^ e ||^^- 
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In order to deal with the boundary integral / gld y gl\ y=0 , one may apply the boundary 

./ 

conditions (4.1) { and (4.3)^ to justify that 



dy9%=o = #V + ^d x U - a"g% =0 . 
This boundary condition allows us to reduce the order of d y g € s , and hence, using the simple 



trace estimate (5.15) and the facts that co e \ y= o > 6 and ||a e ||£oo < 5 , one may prove that 



(5.44) 



g £ s d y gl dx\ 



< z \\(i+vr^gi\\h+c6{i+\\mi> 



.<E II 2 



^+L7||^ + V||i 2(T) . 



Substituting (5.44) into (5.43) we obtain (5.35) 



Proof of (5.36) 



Integrating by parts (cf. remark I5T71) , and using d x u e + d y v e = 0, we have 
(5-45) JJ(1+ y)^g e s {u e d x gl + v%gt} = 7 J J {I + yf^v^gf. 



which and inequality (B.9) imply inequality (5.36) 



Proof of (5.37) 



Since u e e C([0, T]; #*J 4,7 ), it follows from the definition of H s a J a that (1 + y) a u e > 5 
(1 + y y +a2 D a u e \ < 5- 1 for all \a\ < 2. Thus, we have ||(1 + y)d x a e \\ L °° < 5~ 2 + <T 4 

< 5~ 2 , and hence, 



s+4,7 



and 
and 



(5.46) 



2e 2 



(i + y) 2y g e s {d: +1 ^-u) 



d r co ( 



10 



d s x (u e -U)}d x a ( 



< 2e 2 c 5 \\(i + yrgt\M\\(i + vV' 1 ^ 1 ^ - u)\\„ + iki + yy-'d^ - U)\\»}. 

Now, we require the following inequality: 
(5.47) 



Assuming (5.47) for the moment, we can apply it and proposition IB. 31 to (5.46), and obtain 



2e< 



(i + y) 2 XR + V - £0 



M + li 



d s x (u e - U)}d x a ( 



< e 2 C mS {\\u e \\H°<-> + \\d s x U\\ LHJ) + \\d s x +1 U\\ LHJ) + ||(1 + yyd K gl\\v}\\u, e \\ H .« 

1 



which implies (5.37) by Cauchy's inequality and the inequality ||9*Z7||x,3(t) < — \\d x + C7||z: 2 (t)- 



To complete the justification of (5.37) , we have to verify (5.47) as follows. 
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Since 8 < (1 + y) a u e < 5~ x and u e \ y= o = 0, applying part (ii) of lemma iRTl we have 



(1 + 2/)" 



IKi+^-^+V-^IU^ r 1 

(5-48) < , „,,. 

It is worth to note that 



L 2 



+ 1 



^IU 2 (T) + 



L 2 



a, ( ^ ) = g e s+1 = d x gt + d x a e %(u< - U), 



ijj'- 



so by (5.48), we have 

\\(l + yy- l d s x +l (u*-U)\\ L > 
< C w {||^. +1 t/|| L2(T) + ||(l + y)^|| i2 + \\(l + y)d x a e \\ L ^\\(l + y y- 1 d s x (u e -U)\\ L 2} 

which implies inequality (5.47) because ||(1 + y)d x a e \\L°° < S~ 2 + 5~ A . 



Proof of (5.38 



The inequality (5.38) follows from Ho^g^Hloo < 5 +5 and the definition of ||ct/||# s >7. 



Proof of (5.39) and (5.40) 



Both inequalities (5.39) and (5.40) follows from the Holder inequality and proposition IB. 31 



Proof of (5.41) 



For j — 2, 3, • • • , s — 1, using proposition IB.3I and ||(1 + y)a £ \\L°° < 8 2 , we have 

//(i + »)^-v { ^.-».« s ; 



< 



l + y 



{||(i + y y +1 did y ^\\ L2 + 11(1 + y )a*\\ L ~ ||(i + y) WIIWIKi + yVgl 



\L 2 



L°° 



(5.49) 

<C Sl7)CT 4lk1lH|- + ll^^ll^(T)}|k e ||^,T. 

When j = 1, using proposition IB. 31 and ||(1 + y)a e \\L°° < 5 -2 again, we have 

J J (i + y ) 2 ^r VR.<v e - a e a^}^ 

S£- V + yd s x U 



< 



+ 



l + y 



L°°(T) 



L 2 



;i + ^ +1 9^i L oo + ||(i + ^^11^11(1 + ^^1^) 



;i + yy +1 d x d y u e \\ L 2 + || (1 + y)a e |U- 1| (1 + yXWlM } || (1 + y)^|| L2 
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(5.50) 



s,7 . 

9 



Combining estimates (5.49) and (5.50), we prove inequality (5.41) 



Proof of (5.42) 



For any j — 0, 1, • • • , s — 1, by proposition IB.3I and ||(1 + y)a e \\L°° < 5 2 , we have 
jj{l + y)^g^di{u^-U)d^ +1 U 

< ||(1 H-y^^II^IKl -h^^Hi.-IICl -hj/^-^Cx** — COII^II^~ J * +1 ^ r |l-c-cr> 

< C7 SiW ||^ +1 f/|| L ^ (T) {||a;l^„ + ||^f/|| L2(T) }||a;l^, 



which implies inequality (5.42) 



□ 



5.1.3. Weighted H s Estimate on uj e . The aim of this subsubsection is to combine the 



estimates in proposition 15.61 and proposition 15.101 to derive the growth rate control (5.T 
on the weighted H s energy of a/. 

According to propositions 15.61 and I5.l0| we know from the definition of || • Wh'^ that 

d II 6112 
— ||Ct> \\jjs,~1 



dt 



< C S ^, 5 {1 + \\^\\ Hl , + \\dlU\\ L ~ {T) }{\\u:*\\ Hr + \\dl +1 U\\ L ~ m }\\^\\ Hri 

3 

2 

+ C^ 5 \\d s x U\\ |oo( T )||w e ||^,7 + C S;7;CT> 5{1 + ||a; e || H |,7} s 2 ||w e ||^,7 + C s \\9tdxP e \\ 2 H s-2i^ 



<C S + C S , W {1 + ||^ +1 f/|| 4 LO o (T) } + C s Ud x p 



1=0 



e||2 

H s ~ 2l (T) ' 



/=0 



and hence, it follows from the comparison principle of ordinary differential equations that 
\\" e {t)\\h 



rs,7 

1 g 



2 

" s-2 



< [Ikoll^ + J Q F(r)drj jl- (- - 1)C S , W (JM^ + J q F(r)drj I 
as long as the second braces on the right hand side of the above inequality is positive, where 



F : [0, T] 



is defined by (5.2) This proves inequality (5.1) 



5.2. Weighted L°° Estimates on Lower Order Terms. In this subsection we will 
derive uniform (in e) weighted L°° estimates on D a u e for \a\ < 2 by using the classical 
maximum principle. The key idea is to "view" the evolution equation of D a u e as a "linear" 
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parabolic equation with coefficients involving higher order terms of u e ,v e and w e , which can 
be controlled by proposition IB. 31 provided that ||w £ ||h s, t' < +00. 
More precisely, we will prove part (ii) of proposition 15.31 as follows. 

Proof of part (ii) of proposition [3731 This proof, based on a simple application of the classical 
maximum principle for parabolic equations, will be divided into two steps. In the first step, 
we will derive weighted L°° controls on D a u e by using the maximum principle stated in 
appendix [El These controls will rely on the boundary values of D a u e at y = 0, so we will 
also derive estimates on the boundary values of D a uj e by using Sobolev embedding or growth 
rate control argument in the second step. 

Step 1: (Maximum Principle Argument) 

First of all, let us derive a L°° estimate on the I := |(1 + yY +a2 D a uj e \ 2 as follows. 

\a\<2 

For notational convenience, let us denote, for |a| < 2, 

B a := (1 + y) a+a2 D a u e 

which is our concerned quantities. By a direct computation, B a satisfies 

(5.51) {d t + u e d x + v e d y - t 2 8 2 x - 8 2 y }B a = Q a 8 y B a + R a B a + S a , 

where the quantities Q a , R a and S a are given explicitly by 

2(a + a 2 ) cr + a 2 e (a + a 2 )(cr + a 2 + 1) 

:= > K a '■= ~, — ; v + 



{0 if a = (0,0), 

" E (fj + yf 2 D^B a _? +ei + (1 + yf^D^B a ^ +e2 } if \a\ > 1. 

Here, e\ := (1,0) and e 2 := (0,1). Using proposition IB. 31 we have the following pointwise 
controls on Q a ,R a and S a : for |a| < 2, 

" IQol < C a , \R a \ < C mS {\ + \\oo\\ Hr , + ||^.f/|| L2(T) }, 

I'S'al < Cs,7, CT ,5{||w||// g s ^ + ll^x'^IU 2 } z2 i\ B a-l3+ei\ + \B a -p +e2 \}, 

0</3<a 

where C a and C ar/ a - $ are some universal constants which are independent of the solution u e . 
Let us recall from the definition that / := \B a \ 2 , so using (5.51) and (5.52), we have 



(5.52) 



|o|<2 

,2q2 



{d t + u e d x + v £ d y -e'd' x -d' y }I 
- 2 {t 2 \d x B a \ 2 + \d y B a \ 2 } + 2 {QaB a d y B a + R \B a \ + S a B a } 

\a\<2 \a\<2 
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<C S)7jCT)(5 {l + ||a; e || if |,7 + ||^C/|| X 2}J. 
Applying the classical maximum principle for parabolic equations (see lemma IE. II for 



instance) to the quantity J, we have, after using the definition (5.4) of G 
\\I 

(5.53) 



L°°(TxR+) 



re[o,t] 

Next, we are going to derive a lower bound estimate on -B(o,o) := (1 + y) aujt ■ To do this, 
let us recall that -B(o,o) satisfies 

{d t + u e d x + (v £ - Q(ofl))d y - e 2 d 2 x - d 2 }B {m = #(o,o)-B(o,o)- 



Using the classical maximum principle (see lemma |R21 for instance) and (5.52) we obtain 
(5.54) 



min(l + yYuoHt) 

TxR+ 



> (1 - C. W {1 + G(t)}te c ^ s{1+Gmt ) min{ min (1 + y)°u) , min co% =0 }. 

V ' ' TxR+ [0,t]xT U 

Step 2: (Controls on Boundary Values) 



According to inequalities (5.53) and (5.54) , we have already controlled the underlying 



quantities I and -B(o,o) '■= (1 + y) a u) € by their initial and boundary values. However, their 
boundary values are not given in the problem, so we will estimate them in this step. 

\y=o 



In order to control J|„=n := \D a u € \ y= o\ 2 , we will apply Sobolev embedding argument 



M<2 

and growth rate control argument in the cases s > 4 and s > 6 respectively. Combining the 



boundary estimates on / with (5.53) , we will finally obtain inequalities (5.3) and (5.5) 
To derive inequality |(5.3) we first apply lemma lB~2l to obtain 



\I\ y =o\\ L ~ m < Yl H / IU~(TrxR+) < 3C 2 J2{\\D a co e \\ 2 L2 + \\d x D a co e \\ 2 L2 + \\d 2 D a ^\\ 2 L2 } 

\a\<2 \a\<2 



< 6C z \\lo 



2 1 1 . ,e||2 



trig 



which and inequality (5.53) imply inequality (5.3) 



To derive inequality (5.5), let us begin by writing down the evolution equation for 



at y = 0: for any \a\ < 2, 
(5.55) d t D a uj 
where the term E a is given explicitly by 





e 2 d 2 x + d 2 )D a u e + E a 



y=o 



if «2 = 

if ai = and a 2 > 1 
d x u e \ 2 - u e d 2 x uj e if a = (1,1). 
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Here, the derivation of (5.55) is just a direct differentiation on the vorticity equation (4.3) 
as well as using the boundary condition u e \ y= o = v e \ y= o 
IR31 if s > 4, then 



0. Furthermore, by proposition 



(5.56) ||-E'«|?/=o||l° o (t) < Cs,7l|a/||/^'T- 

In addition, if s > \a\ + 4, then by proposition IB .31 again, for e G [0, 1], 

(5.57) \\(e 2 d 2 x + d 2 v )D a u% =0 \\ L ~ m < C s J^\\ H ,-r. 



Therefore, using (5.55) - (5.57) and inequality (B.10), we have, for s > 6, 

l|dt/L=olU» (T) < CsA 1 + IMk^liK" 2 



and hence, by direction integration and definition (5.4) of f2, we obtain 

(5.58) ||/(t)|, = olU-(T) < \\i(o)\ y =o\\L~(T) + c s , 7 {i + n(t)}n(t)H. 



Combining estimates (5.53) and (5.58), we prove our desired estimate (5.5) 



Lastly, it remains to show inequality (5.6) Let us begin by deriving an estimate on u e 



l^o- 
Using identity (5.55), inequality (5.57) and definition (5.4) of fl, we have, for any s > 4, 

||^o; <e | ?/= o||loo(t) < C s ,7^(£), 

so a direct integration yields 

(5.59) minc/^l^o > minwo^o — C s ^Q(t)t. 



Combining estimates (5.54) and (5.59), we prove (5.6) 



□ 



6. Proof of the Main Theorem 



The purpose of this section is to complete the proof of our main theorem 12.21 In other 



words, we will prove existence and uniqueness to the Prandtl equations (2.1) in the subsection 
16.11 and 16.21 respectively. 

6.1. Existence for the Prandtl Equations. In this subsection we will construct the 



solution to the Prandtl equations (2.1) by passing to the limit e goes to + in the regularized 



Prandtl equations (4.1) Our proof will be based on the uniform (in e) weighted estimates 
derived in proposition [531 Using these estimates, we will derive uniform bounds and lifespan 
on a/, and then prove convergence of u £ and consistency of the limit u as follows. 

Uniform Bounds and Life-span on u e 



According to proposition 15.11 a solution uj € to the regularized Prandtl equations (4.1) 
exists up to a time interval [0, T S)7)O - )( 5 i(EiW0i t/], which may depend on e as well. However, in 
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proposition 15.31 we have already derived uniform (in e) estimates on u e , so one may apply 



the standard continuous induction argument to further solve the Prandtl equations (2.1) up 
to a time interval which is independent of e. As a result, we have 



Proposition 6.1 (Uniform Life-span and Estimates for u € ). In addition to the hypotheses 
of proposition \5.1\ when s = 4, we further assume that S > is chosen small enough 



such that the initial hypothesis (2.7) holds. Then there exists a uniform life-span T 



T(s, 7, a, 5, || c^o 1 1 U) > 0, which is independent of e, such that the regularized vorticity 
system \(4~3j\ - [(44)1 has a solution cu e G C([0, T];H*'J) D C fl ([0, T]; H s ~ 2,1 ) with the following 
uniform (in e) estimates: 

(i) (Uniform Weighted H s Estimate) For any e G [0, 1] and any t G [0, T], 
(6-1) \\^(t)\\ H s. <4\\u \\ H s„. 

(ii) (Uniform Weighted L°° Bound) For any e G [0, 1] and t G [0,T], 



(6.2) 



\(i+yy +a2 D a u;*(t)\ 

\a\<2 



1 

< — . 

- 5 2 



L°°(TxR+) 

(iii) (Uniform Weighted L°° Lower Bound) For any e G [0, 1] and t G [0,T], 



(6.3) 



mm(l + y) a u € (t) > 5. 

TxR+ 



Proof. The uniform life-span T := T(s, 7, a, 5, \\uo\\h s -/, U) can be guaranteed by the uniform 
estimates (6.1) - (6.3), so it suffices to justify them. Indeed, the life-span T can be taken as 



min{T l5 T 2 , T 3 } where T 1; T 2 and T 3 will be defined below. 



(i) According to the definition (5.2) of F and the regularized Bernoulli's law (4.2) 

-+1 

II-^IIl 00 < Cs,7,cr,5{l + II^^IIh s - 2 '+ 2 (T)} 2 — Cs^^sMu 



1=0 



f +1 c» II, , 1 1 2 

where M v := sup{l + > ||^[/||^ s „ 2 (+2 m } 2 < +00, so if we take T t := min < — f— , 



1-2 



-s+2 



1=0 



-2 ' 



7,cr,(5 1 1^011/^.7 



then by inequality (5.1) 



estimate 



(6.1) holds for all t G [0,Ti]. 



(ii) When s > 6, using part (i) of proposition 16.11 we know from the definition (5.4) of 
fl and G that for any i 6 [0, Ti] where Ti is defined in part (i), 



(6.4) 



M(t) < 4\\u) \\ H s 



and G(t) < 4||w ||h«.'» + M u =■ K. 
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Thus, if we take T 2 := min < 7\, 



In 2 



■, then 



645 2 C\ 7 (1 + M\u \\ H sn) \\cj \\ 2 Hl n ' C m5 (l + K) 
using inequality 1(5.5)1 and the initial assumption > |(1 + y) aJra2 D a coo\ 2 < — -, we have 



\a\<2 



the upper bound (6.2) for all t G [0,T 2 ]. 

estimate 



When s = 4, using inequality (5.3) 



(6.4) and the initial hypothesis (2.7) , we also 



have the upper bound (6.2) for all t G [0,T 2 ]. 



(iii) Let us take T-t 



min < Ti, 



In 2 



8C S 



Then 



using inequalities (5.6) and (6.4) , we know that the lower bound (6.3) holds for all t G [0, Tj 



□ 



Convergence and Consistency 

Using almost equivalence relation (A.l) and uniform weighted H s estimate (6.1), we have 

(6.5) sup (||w e ||ir.,i + \\u e - U\\ H .,i-i) < C,, 7 ,« r ,j{4||wb||fl;^ + II^^IU 2 } < +°°- 

0<t<T 



Furthermore, using evolution equations (5.29) , uniform H s bound (6.5) and proposition IB. 31 
one also find that dtu e and dt(u e — U) are uniformly (in e) bounded in L°°([0, T]; H s ~ 2,1 ) and 
L°°([0, T]; H s ~ 2,1 ~ ) respectively. By the Lions-Aubin Lemma and the compact embedding 
of H sn in Hi oc stated in lemma |6~2| we have, taking a subsequence if necessary, as — > + , 

u ek ± u in L°° ( [0, T] ; H s ' r ) and w £fc ->win C([0,T]; HQ , 

u tk - U ^ u - U in L°° ( [0, T] ; iP' 7 " 1 ) and u efc -»■ u in C( [0, T] ; HQ , 

for all s' < s, where w G L°°([0, T]; fP' 7 ) n p| C([0, T];Hf 0C ), u-U G L°°([0, T]; IP' 7-1 ) n 

s'<s 

C([0,T]; #f oc ) and w = (9 y M. Using the local uniform convergence of d x u 6k , we also have 



(6.6) 



s'<s 



the pointwise convergence of v tk : as e& — > 



(6.7) 



- / d x u ek dy -> 



c^w =: 



Combining (6.6) - (6.7), one may justify the pointwise convergences of all terms in the 



regularized Prandtl equations (4.1) L — (4.1)^ Thus, passing to the limit — » in (4.1' 



(4.1) j and the regularized Bernoulli's law (4.2), we know that the limit (u,v) solves the 



Prandtl equations (2.1) L — (2.1) l with the Bernoulli's law (2.3) in the classical sense. 

Lastly, in order to complete the proof of consistency, it remains to justify that u G 
L°°([0, T]; H^J) and the matching condition (2.1).. Since D a u 6k converges to D a u 
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(6.8) 



J2\(i + y y +a w a u^\ 2 <^ 



\a\<2 



and 



;i + y) ff w e * >6, 



we deduce that (6.8) still holds for u, and hence, u G L°°(0,T; if 
Also, by the Lebesgue's dominated convergence theorem, 



bj dy 



lim 



tu efc dy = U 



which is equivalent to the matching condition (2.1) ; because ui = d y u > 0. 
To complete the proof of existence, let us state and prove the following 

Lemma 6.2. Let s be a positive integer, 7' > and M < +00. Assume 

(6-9) ||r|Uy <M 

for all e G (0, 1]. Then there exist a function f G H sri and a sequence {ek]ken Q (0, 1] with 
lim 6^ = + such that as e& — > + , 

k— >+oo 

(6.10) f k H -?' f and f ek H -k f for all s' < s. 

Proof of lemma \6.B . First of all, let us mention that H s ' y has an inner product structure: 

1 Jo Jt 

\a\<s 



so the uniform bound (6.9) implies the weak convergence of f ek in (6.10) via the Banach- 
Alaoglu theorem. 

Next, by the definition of || • H^.y, ||/ e ||// s < II/^IIh 8 ^' — M. This implies the local H s 
norm convergence in 



(6.10) for all s' < s because of the standard compactness of H S (T x ' 



□ 



Finally, let us end this subsection by giving the following 



Remark 6.3 (Life-span for U = constant). In the special case that U > is a constant, one 
may prove that the life-span T stated in our main theorem !2.2l is independent of the constant 
value of U, that is, T := T(s, 7, a, 5, ||wo||h s ^)- The reasoning is as follows: 

When U = constant, it follows from the regularized Bernoulli's law (4.2) that d x p e = 0, 
so by definitions (5.2) and (5.4), we have F = C S7(T 5 and G(t) = Q(t) = sup Uu^H^.-r where 

[0,t] 



all of F, G and f2 are independent of U. As a result, all of our weighted estimates (5.1) 



(5.3) , (5.5) and (5.6) are independent of U. Therefore, one may slightly modify the proof of 



proposition 16. II to show that uniform weighted estimates (6.1) - (6.3) hold in a time interval 
[0, T St y iat s t \\u) \\ Hsn ), which is independent of U. According to our proof of convergence and 
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consistency in subsection 16. 1\ we can solve the solution (u,v) of the Prandtl equations (2.1) 
in the same time interval, and hence, the life-span T stated in the main theorem 12.21 is also 
independent of U. 

6.2. Uniqueness for the Prandtl Equations. The aim of this section is to prove the 
uniqueness of H^J solutions constructed in subsection 16.11 To show the uniqueness, we will 
generalize the nonlinear cancelation applied in subsubsection 15.1.21 to the L 2 comparison of 
two H^J solutions. This motivates us to consider the quantity g below. 



Specifically, the uniqueness of H^'] solutions to the Prandtl equations (2.1) is a direct 
consequence of the following L 2 comparison principle. 



Proposition 6.4 (L 2 Comparison Principle). For any s > 4, 7 > 1, a > 7 + - and 5 £ (0, 1), 



let (ui,Vi) solve the Prandtl equations (2.1) with the vorticity Ui := d y u,i £ C([0, T]; H s jJ) fl 



C 1 ([0,T];iJ s 2 ' 7 ) fori = 1,2. Define g := U\ — U2 + ® yU>2 {u x — u 2 ) ■ Then we have 



ijj 2 



(6.11) 



2 

L 2 



~ 1 1 2 



\d y g\\ 2 L 2 < ||^(0)||i 2 +C 1>aA u,u I 11:711^ 



where the constant C lja} s jU) ,u depends on 7, a, S, HcjiH^^, H^H^^ and ||c^^||l 2 (t) only. 



Applying the Gronwall's lemma to (6.11), we obtain 

wmwb < \\m iii 2 e c — * 

which implies g = provided that ui\ t =o = U2\t=o- Since U2d y 
(6.12) ui-u 2 = qu) 2 



Ui — u 2 
0J2 



g = 0, we have 



for some function q := q(t,x). Using the Oleinik's monotonicity assumption U2 > and the 
Dirichlet boundary condition u,j\ y= o = for i = 1,2, we know via (6.12) that q = 0, and 



hence, u\ = U2- Since Vi can be uniquely determined by Ui, we also have v\ = V2- This 
proves the uniqueness of H s a J solutions. 

In the rest of this subsection, we will prove proposition 16.41 as follows. 



Proof of proposition \6-4\ Let us denote (u, v 

u \ 



[Ui,Vi, 



(u 2 ,V2) and 02 '■= ^ yU2 . Then one 



UJ 2 



may check that g = oj — a 2 u = u 2 d y ( — I and satisfies 



(6.13) 



(d t + u\d x + vidy - d y )g = -2ujd y a 2 - u{ud x a 2 + vd y a 2 + 2a 2 d y a 2 }. 



To derive the L 2 estimates on g, let us first recall that we define the cutoff function 



y 



Xn(y) '■= x(~e>) f° r an y -R > 1) where x ^ C™([0, +00)) satisfies the properties (4.7) Then 
R 
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Xr has the following pointwise properties: as R —¥ +00, 

\Xr\<^0 + and \ Xr \ < O(^) -+ 0" 



For any t G (0, T], multiplying equation (6.13) by 2xng, and then integrating over [0,t] x 
T x IR + , we obtain, via integration by parts, 

XR9 2 {t) dydx - J J XR9 2 \t=o dydx 

= -2 / // XR\d y g\ 2 -2 / / gd y g\y = o dx - 4 / // XRgud y a 2 
Jo J J Jo Jt Jo J J 

— 2 J J J XR9u{ud x a 2 + vd y a 2 + 2a 2 d y a 2 } + 7£i + 7£ 2 

< - 2 / // x it |<9^| 2 -2 / / gd y g\ y=0 dx + 4||<9,,a 2 || LO o / I^IUHI^IU 2 
Jo J J Jo Jt Jo 



(6.14) +2||(1 + ?/){£i<9 x .a 2 + vd y a 2 + 2a 2 d y a 2 }\\ L ^ / 
where the remainder terms IZi are defined by 



a 



1 + y 



+ ^i+^ 2 , 



L-' 



Ki:= j jj ^ vi ~ g2 and U2 - = So SS x *~ g2 



O^JJL 



Now, the first technical problem is to deal with the boundary integral 

Since d y g\y=o = — a 2 g\ y= o, we have, after applying the simple trace estimate (5.15) 
(6.15) 



gd y g\ y =o dx. 



gd y g\ y=0 dx 



Jt 



< 



7, J jj XR\d y g\ 2 + C{||a 2 || LO o + Hoallioo + ||9 v a 2 ||r<»} J J J <f- 



Furthermore, since u 2 G H^'], it follows from the weighted L°° bounds on u 2 that 



||(l + y)a 2 || LO c < cT 2 
1 + y)d x a 2 \\ L ™, ||(1 + y) 2 d y a 2 \\ L oo < 5' 2 + (T 4 , 



(6.16) 

so by proposition IB.3| 

(6.17) \\(l+y){ud x a 2 +vd y a 2 + 2a 2 d y a 2 }\\ L o C < l7 7i<Tj(5 {1 + ||o;i|| h 4, 7 + ||w 2 || H 4, 7 + ||<9 4 £/|| L 2 (T) }. 



Substituting (6.15) - (6.17) into (6.14), we obtain 



XRg (t) dydx - XRg \t=o dydx 



(6.18) ^-JJJ **\ d v9\* + C * I 



^7,<j,i5, ||oji I 



„4,7.II W 2|L4, 7 ,||9|C/|| L 2 
n g "g 



g\\L* s IMU 2 



l + y 



L 2 



+ n 1 + n 2 . 
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Next, we emphasize that both uj and 
Claim 6.5. 

(6.19) \\Cj\\ L 2, 



u 



1 + y 



u 



can be controlled by g, namely, 



L 2 



The proof of claim 16.51 is very similar to that of lemma \A.2\ so we will only outline it at 



the end of this subsection. Assuming claim 1631 for the moment, we can substitute (6.19) 
into 



(6.18) to obtain 

XR9 2 (t) dydx 



(6.20) 



Xr9%=o dydx 



< 



Xn\dy9\ 2 + C-y,a,6,\\ 



w lH„4, 7 ,||^2||„4, 7 ,||C'JU|| i 2 

9 9 JO 



r ,\mu\\. 



\g\\ 2 L 2+Ki + K 2 . 



Finally, both integrands of 7Zi and 1Z 2 can be controlled by a multiple of g 2 , which belongs 
to L\[0,T];T x R + ), so applying Lebesgue's dominated convergence theorem, we have 



(6.21) 



lim IZi 

R-^+oo 







for z = l,2. 



Using monotone convergence theorem and (6.21) , we can pass to the limit R — > +oo in 



(6.20) to obtain (6.11 



5 



Lastly, we will justify claim IBT51 as follows. 
<, \6.5l Using triangle inequ 
, so it suffices to control u. 



Proof of claim WH Using triangle inequality and (6.16) L , we have ||cj||l2 < \\g\\L 2 + 
u 



-2 



1 + y 



L 2 



Since 5 < (1 + yf^i ^ ^ 1 an d ^1^=0 = 0, applying part (ii) of lemma [B?T| we obtain 



u 




1 + y 


L 2 



[i + y) 



-CT-1 



Lu'2 



<c a , ( 



L 2 



u 

^2 



< Ca,<5||fi , ||L 2 



L 2 



because g = uizdy 



u 



0J2 



□ 
□ 



7. Existence for the Regularized Prandtl Equations 



The aim of this section is to solve the regularized Prandtl equations (4.1) , or equivalently 



its vorticity system (4.3) - (4.4) In other words, we will prove proposition 15.11 according to 
the plan described in section HI However, we will only sketch our proof because the method 
for solving intermediate approximate systems (4.3) - (4.4), (4.5) - (4.6) and (4.8) - (4.9) is 



standard. Before we proceed, it should be also remarked that we will solve the approximate 



systems (4.8) - (4.9) , (4.5) - (4.6) and (4.3) - (4.4) with a decreasing order of regularities 



The main reason of this technical arrangement is to derive our estimates in a rigorous way so 



EXISTENCE AND UNIQUENESS TO THE PRANDTL EQUATIONS 



35 



that we can differentiate the intermediate equations pointwisely and have enough pointwise 
decay at y = +00 according to prop IC.ll 



7.1. Solvability of Inhomogenous Heat Equation. In this first subsection we will solve 
an inhomogenous heat equation in the weighted space H s a 'J. This existence result will be 



applied to solve the linearized, truncated and regularized vorticity system (4.8) - (4.9) in the 
next subsection. 

Let us consider the following inhomogenous heat equation: for any e > 0, 
d t W + F R = e 2 d 2 x W + d 2 W in [0, T] x T x R + 
W\ t=0 = W on T x R + 

dyW\ y=0 = d x p € on[0,T]xT 
where W is an unknown, Wq and d x p £ are given initial and boundary data, Fr is a given 



(7.1) 



inhomogenous term with compact support in [0, T] x T x [0, 2R\. Since (7.1) is just a standard 



inhomogenous heat equation, we can solve it by classical methods and obtain the following 

Proposition 7.1 (Existence of Inhomogenous Heat Equation). Let s > 4 be an even integer, 
7>l,<r>7+i <fe (0,i) andee (0,1]. //W G 7 and suppi^C [0, T] x T x [0, 2i?] ; 

then there exist a time T := T(s, 7, a, S, R, ||Wo||fl-»+8,7, H-FrIIc 2 ) > and a solution W G 
L7([0,T];#*J 8 ' 7 ) n L7°°((0,T] x T x R + ) to the inhomogenous heat equation |(7l) . 

Furthermore, we have the following pointwise decay aty = +00; for any I — 0, 1, • • • , — +4, 
for any \a\ < s — 21 + 9, 



(7.2) 



d[D a W 



o((i + y y°- ° 2 ) 

IT + (2 |a| + 2i-2_ 1)7 

0((l + y) 



if\a\ + 2l<2 
if2<\a\ + 2l<s + 9, 



and energy estimate: 

d 
If 



(7.3) 



|W||ls + 8,7 + e2 HI^III' + 8,7 +IH^I 



12 

ls+8,7 



< c 



s,7 I 



ls+8,7 



+ c s \\\w\ 



\s+8,-y\ 



\Fr\ ||s+8,7 



|2 

ls+8; 



where the norms 



and 



are defined in definition B.4 



Outline of the Proof. Using the method of reflection and Duhamel's principle, one may 
express the unique global-in-time C([0,T] x T x M+) n C°°((0,T] x T x R + ) solution to 



(7.1) by an explicit solution formula 

(7.4) W = K Wo + K dxpt + K Fr 

where the terms K Wo ,Kg xP ^ and K Fr can be written explicitly by using the Gaussian (i.e., 
the heat kernel), and depend on Wo,d x p e and Fr respectively. Since the solution formula 
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(7.4) is explicit, based on the properties of the Gaussian, one may prove the following two 
facts: as y — > +00, 

(i) both Kq xP * and Kp R decay exponentially fast; 

(ii) the term D a K Wo < (or >)(1 + y)~ K provided that D a W does. 



Using quantitative versions of facts (i) and (ii), one can justify by (7.4) that W fulfills all 
weighted L°° controls for H^ 8,1 within a short time interval [0 } T sn ^s,R,\\Wo\\H s rf,\\F R \\ c2 ]- 

Furthermore, applying proposition IC.ll with s' = s + 8, we know that Wq satisfies the 
pointwise decay (7.2) for / = and |a| < s + 9, and hence, W does. Using the heat equation 



(7.1) L repeatedly, we also obtain the pointwise decay (7.2) in our desired ranges of / and a. 



Finally, it remains to show the energy estimate (7.3) but its proof just follows from 



standard energy methods, so we will omit the proof here. However, during the estimation, 
one requires to apply an integration by parts in the ^/-direction to deal with the operator dy, 
so we would like to give the following two remarks on the boundary values of W: 

(I) (Boundary Values at y — 0) The boundary values of W as well as its derivatives can 
be reconstructed by using the boundary reduction formula 

jfe-i 



d 2k+l W 



y=0 = (d t - e 2 dl) k d x f + - e 2 d 2 x ) k -^ +1 F R \ y=0 , 

j=0 

which reduces the order of the boundary terms so that we can control the boundary 
integral at y = via the simple trace estimate (5.15) 
(II) (Boundary Values at y — +00) All boundary terms of W as well as its derivatives 

+00 



required for deriving energy estimate (7.3) actually vanish fast enough at y 



because of the pointwise decay estimate (7.2) Thus, all required boundary integrals 
at y = +00 are zero. 

□ 



7.2. Solvability of Linearized, Truncated and Regularized Vorticity System. Based 



on the H^l solutions to the inhomogeneous heat equation (7.1) derived in subsection 17. 11 we 



will construct a sequence of solutions to the linearized, truncated and regularized vorticity 



system (4.8) - (4.9) with uniform bounds in this subsection. This sequence of solutions as 
well as their uniform bounds will be the foundation for solving the truncated and regularized 



vorticity system (4.5) - (4.6) in the next subsection. 

Let us begin by defining an iterative sequence {(u n , v n , o; n )} ne p 

(i) u°(t,x,y) := u (x,y); 
(ii) 



as follows: 



(u n ,v n ) 



is defined by formulae (4.9) for all n e N 



111 0J 



n+l is defined to be the C([0, T]; H & a + * a ) n C°°((0,T] x T x 



solution to the 



linearized, truncated and regularized vorticity system (4.8) - (4.9) for all n e iV 
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The natural question is whether the iterative sequence {(u n ,v n ,u n )} n€ ^ is well-defined, 
and the answer is affirmative because of the following 



Proposition 7.2 (Existence of Linearized, Truncated and Regularized Vorticity System). 

1 1 

Let s > 4 be an even integer, 7 > 1, cr > 7+ G (0, -),e G (0,1] and R > 1. // 



ujq G -f^ttf 2 ' 7 an d SU P II 1^1 lls+9,00 < +00 where 
t 

there exist a uniform life-span T := T(s,j,a,5,e,x, R, \\^o\\h b + 8 ^, su P II \U\ || s +9,oo) > which 

t 

is independent of n, and a sequence of solutions {u n } n &i Q C([0, T]; H^ 8 '" 1 ) fl C°°((0,T] x 
T x I 



is defined in definition B.4, then 



to the linearized, truncated and regularized vorticity system (4.8) - (4.9) 



Furthermore, the pointwise decay estimate (7.2) holds for W := u n for all n G N and we 



have the following uniform (in n) energy estimate: for all n G N and for all t G [0, T], 



(7.5) 



n\\\2 



I s+8,7 



t pt 

s+8i7 + / \\\d y UJ 



\\d x u 



nl ||2 



n| ||2 

s+8,7 



< Qs+io(II w oIIh s + 8 , 7 ) 



where the norm 



||| s ' , 7 is defined in definition B.4 and Qi is a degree I polynomial with 
non-negative coefficients which depends on C Sj7iX and |||C/||| s+8)00 only. 



Outline of the proof. For a given u n G C([0, T]; H s a + 8n ) n l7°°((0, T]xTxl + ), the local in 
time solvability of u n+1 in the same function space follows directly by applying proposition 
Owith W := u n+1 and F R := xn{u n d x Lo n + v n d y u n }, although the life-span T may depend 



on n a priori. However, the unform (in n) energy estimate (7.5) guarantee the unform (in 



n) life-span T by the standard continuous induction argument, so it suffices to prove (7.5) 

In order to derive the energy estimate (7.5) , we have to control || \Fr\ Hs+7,7 and || \F R \ \\ s+8 >7 
for Fr := x/j{M n 9 x a; n + v n d y uj n }. Using the triangle inequality, proposition IB.5I and 
proposition IB.6| one may check that 

|-Fr|||s+7,7 - ^'s,7,xlll a;n |lls+8, 7 + 111^1 



(7.6) 



|2 

I s+8,00 



ill || s+8,7 



where the norms 



— Cs,7,x,m{|| l 1 ^™! lls+8,7 + || |^| ||s+9,oo} ' { II I^^IH s+8,7 + II I^^M II s+8,7} 

and III ■ I ||s' 00 are defined in definition IB. 41 



| s',00,7 



Applying inequalities (7.3) and (7.6) , one can easily show that as long as || \u n \ 



|s+8,7lt=0 



< 



L, there exists a uniform (in n) time interval [0,T a)7 tT ^ e)X( jj. 



su Pt lllf HU+9,c 



(7.7) 



\U) 



nl 1 1 2 



I s+8,7 



\\d x u 



nl 1 1 2 



I s+8,7 



+ 



\\dyU 



n I II 2 < AT 2 
s+8,7 - ^ L 



A such that 



for all n G N 



because both (7.3) and (7.6) are independent of n. Therefore, it remains to derive a uniform 
(in n) control on the initial data |||a; n ||| s -|-8 i7 |t = o. 



To estimate \\\u r 



(71 



| s+8j7 | t=0; let us first state without proof the following fact: 
l|d t w n ||# s -2!+8, 7 | t= o < ?Vi(ll w o||ip+ 8 ^) 
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where Vi + i is a degree I + 1 polynomial defined by 

2-1 

V X {Z) := Z and V l+1 := V t + C s ^ x ^(Vj+i 

j=0 



WW 



s+8,oo)'Pl-j 



for all I > 1. 



The fact (7.8) can be proved by induction on (n, I) together with the following estimate: 

i 

(7.9) Y n+ 1,1+1 < Yn+l,l + Cs,7,x / J (^n,j + || |^| ||s+8,oo)^n,/-j 

j=0 

where Y^j := ||d t u; n ||# s -2i+8, T |t = o. The derivation of (7.9), which is based on (4.8) L , 
proposition IB .51 and proposition IB.6| will be left for the reader. 

Combining estimates (7.7) and (7.8), we show the uniform energy estimate (7.5) for 

5+4 



Q s+10 := 4j2 V li- 



1=0 



□ 



7.3. Solvability of Truncated and Regularized Vorticity System. The aim of this 



subsection is to construct a solution to the truncated and regularized vorticity system (4.5) 



(4.6) by passing to the limit in its linearized version (4.8) - (4.9), which was solved with 



uniform bounds in subsection 17.21 

In other words, we will prove the following 



Proposition 7.3 (Existence of Truncated and Regularized Vorticity System). Let s > 4 be 



an even integer, 7 > 1, a > 7 + -, 5 E (0, -), e 6 (0, 1] and R > 1. If Uq E H S G 



s+12,7 

25 



and 



sup |||[/||| s+ 9 >00 < +00 where 
t 



S'OD 



is defined in definition B.^, then there exist a time T :- 



T(s,7,cr, 5, e, x,R, ||wo||h s +8.t, sup || |f7||| s+9i00 ) > and a solution lo r E C([0, T]',H*~ 



n 



f]C l ([0,T];H 



to the truncated and regularized vorticity system (4.5) - (4.6) 



1=1 



Furthermore, we have the following uniform (in R) weighted energy estimate: 

(7 10) ^IH W «llls+4,7 + e2 |II^W/j|||s+4,7 + III^-rIIIs+4,7 

< c s>7ie , x {i + 1| M |U >7 + 11 \u\ |U} 2 1| \u> R \ || 2 +4 , 7 + c s {i + \\\u\ \\ s+e } 2 \\\u\ || 2 +6 

and the following weighted L°° estimates: 

(7.11) ||/ii(t)|U«-(TxR+) < (11^(0) |U-(Txr+) + C sn A(t) sup |M|zr.-K,-,t)e c * A W 

[0,t] 



(7.12) min (1 + yYu R (t) > (1 - A(t)te A(t)t )( min (1 + y) a u - C a A(t) sup || 

TxR+ TxR+ [ 0jt ] 



where the norms 
A are defined by 
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I s i7 



are defined in definition B.^, and the quantities Ir and 



I R {t) :=Y,\^ + yY +a2 D a oo R {t)t 

\a\<2 



and 



A(t) := 1 + sup 11^11^+4,7 + sup ||[/|| C 3 (T ). 

[0,t] [0,t] 



Outline of the proof. According to proposition 17.21 the sequence of solutions {u) n } ne $s to the 



linearized, truncated and regularized vorticity system (4.8) - (4.9) has a uniform (in n 
life -span [0,T s ^ tU ^,e,x,R,^>o,u}j i n which || |o; n | || s +s ,7 is uniformly bounded by estimate |(7.5) 
Based on this uniform bound, one may apply the standard energy methods to u n+1 - u n to 

prove that the approximate sequence {co n } n€ ^ is indeed Cauchy in the norm sup ||| • ||| s +6,7 

te[o,T] 

where the time T := T(s, 7, a, 5, e, x, R, \\wo\\s+8,y, SU P II 1^1 lls+9,00) > is independent of n. 



As a result, we can pass to the limit n goes to +00 in (4.8) - (4.9) to obtain a solution 
oJr :— lim bf 1 to the truncated and regularized vorticity system (4.5) - (4.6) Moreover, 

#+4 



cor belongs to C([0,T]; tfj 8 ' 7 ) n f] C l ([0,T]; H s - 2l+8 ^) because co n does. 

1=1 

The uniform energy estimate |(7. 10) follows from the standard energy methods, so its proof 



will be omitted here. It is noteworthy to mention that unlike the estimates in subsection 



17.21 all constants in (7.10) are independent of R. This improvement is based on applying 



integration by parts appropriately to the integral involving the convection term Xr{ u r9x^r+ 



VRdyUji}, but it does not exist in (4.8) - (4.9) because the linearization destroys this structure. 



The weighted L°° controls (7.11) and (7.12) can be derived by the classical maximum 



principle (see lemmas IE. 1 1 and IE. 21 for instance) as in subsection 15.21 We leave this for the 
interested reader. 

□ 

7.4. Solvability of Regularized Vorticity System and Regularized Prandtl Equa- 
tions. In this subsection we will construct a solution u e to the regularized vorticity system 



(4.3) - (4.4) by passing to the limit in its truncated version (4.5) - (4.6) , whose local-in-time 



solvability and uniform bounds are shown in subsection 17.31 Furthermore, we will also justify 



that a/ solves the regularized Prandtl equations (4.1) 



More precisely, we will complete the proof of proposition 15.11 as follows. 



Outline of the proof of proposition [3771 To solve the regularized vorticity system (4.3) - (4.4) 



we first pick any function x with the properties (4.7) in the truncated and regularized 
vorticity system 



(4.5) - (4.6) Then by proposition 17.31 we have a local-in-time solution 



cor to (4.5) - (4.6) and uniform bounds (7.10) - (7.12) on lor. Since the estimates (7.10) 



(7.12) are independent of R, one can show that there exists a uniform time T 



T(s,>y,<T,8,e,\\uo\\H+*.-r,U) > such that {u R } R > 1 C C([0, T]; H^' 1 ) n C 1 ([0, T]; H S+2 ' J ) 
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and ||wfl||ij«+4,7 < Cs,-r,a,s,€,\\ui \\ H3+4 ,^,u f° r an R>1- Therefore, by the standard compactness 
argument, there exist a function u e G C([0, T}; H^' 1 ) nC 1 ([0, T]; H s+2n ) and a subsequence 
{-Rfc}fc G N with lim R k = +00 such that converges to u € in C([0, T]; H£* 2 ) as Rk — >■ +00. 

fc— >+oo 



As a result, we can pass to the limit f4 — » +00 in (4.5) - (4.6) , and prove that u e solves the 



regularized vorticity system (4.3) - (4.4) in a classical sense. 



Finally, we will justify that (u t ,v t ) defined by (4.4) satisfies the regularized Prandtl 



equations (4.1) as follows. 



First of all, the matching condition (4.1) ; , the Dirichlet boundary condition v e \ y= o and 



the initial condition (4.1)^ follows immediately from the formulae (4.4) , loq := d y Uo and 



the compatibility condition (2.2) Then by direct differentiations on (4.4) , we also have the 



incompressibility condition (4.1) 2 and a/ = d y u 
To justify equation (4.1) L , we substitute u e = 



(7.13) 



d y u e into (4.3) L and obtain, via using (4.1) 

d y {d t u e + u e d x u e + v £ uj e } = d y {e 2 d 2 x u e + d y u e }. 



Then one may derive (4.1) L by integrating (7.13) with respect to y over [y, +00) and using 
the following pointwise convergences: as y — > +00, 



(7.14) 



v € u e , d y u e -)■ 

,2q2„,e 



d t u e + u e d x u e - e 2 d 2 u e -»■ d t U + Ud x U - e 2 d 2 U =: -d x p 



The pointwise convergences (7.14) can be shown easily as long as uj e G C([0, T\]H S C 



n 



C 1 ([0, T]; H s+2,y ), so we leave this to the interested reader. 

Lastly, it remains to show the Dirichlet boundary condition u e 



0. To prove this, 

we evaluate the evolution equation (4.1) L at y — 0, and apply the boundary conditions 
v e \ y =o = and (4.3) i to obtain that u e \ y=0 satisfies the viscous Burger's equation: 

e 2 dl(u e \ y=0 ). 



(7.15) 



d t (u e \ y=0 ) + (u e \ y=0 )d x (u e \ y=0 ) 



It follows from the classical uniqueness result for the viscous Burger's equation (7.15) that 
u e \ y= Q = since it does initially according to the compatibility condition (2.2) 



□ 



Appendix A. Almost Equivalence of Weighted Norms 



The purpose of this appendix is to justify the almost equivalence relation (3.2) In other 
words, we will prove the following 

Proposition A.l (Almost Equivalence of Weighted H s Norms). Let s > 4 be an integer, 

7 > 1, a > 7+2 an d ^ e (0> -0- F° r an y u e HajsC^ x ^ + )» we have the following inequality: 
there exist constants c$ and C s ~ a & > such that 



(A.l) 



\u ~ U\\ H .,f-i < C a „ fl r i s{\\u\\ H ',n + \\d x U\\ L 2} 
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provided that uj = d y u, u\ 
and || ■ Wh 3 ' 1 are defined by |(2.5)| and |(3.1 j| respectively. 



y=0 = and lim u = U , where the weighted H s norms || • ||/po 

y— ¥+00 



Proof. Without loss of generality, we only need to prove inequality (A.l) for the smooth ui 
because the full version can be recovered by the standard density argument. First of all, it 
follows from the definition of ||w||fl-«.-r and \\u — U\\h*<i-i that 



Co' 



(A.2) 



k=0 



\ H .n + J2 IK 1 + VV'^liu - U)\\ L 2 < \\u,\\ H .„ + \\U - u\\ H 

< 2{m\h*i + £ IK 1 + yy- ld *( u - u )\\v}. 



s ,7 — 1 



k=0 



Furthermore, applying Wirtinger's inequality in the variable x repeatedly and part (i) of 



lemma |KT| we have 
( s 



EiKi+vr^-toiu^i 



1 + 7T 



-2s 



i + yy-'d^u-U)]] 



k=l 



71 



-2 



L 2 



:i + y y-\u-u)\\ L 2 < -A-iKi + ^n^ < _A_| K .|| i7 . . 



27 - 1 



27 - 1 1 



and hence, there exists a constant G Sil > such that 



(A.3) 



Co' 



\ H .„ + ||(i + y y- l d s x {u - u)\y < \\ U \\ H „ + J2 ll(i + vV^dlin - u)\\„ 



k=0 



< C a J\\u\\ H .n + ||(1 + yy- l d S x {u - U)\\»} 



Therefore, according to inequalities (A.2) and (A.3), it suffices to prove 

(A.4) c 5 \\u\\ Hr < \\u\\ H .„ + ||(1 + V y- l d s x (u - U)\\» < C w {|MUj" + \\d s x U\\ L 2} 

for some constants c$ and C 7i0 -,5 > 0. 

The key idea of proving (A.4) is the following 

Lemma A.2 (L 2 Comparison of d x (u — U), d x u and gt)- Let s > 4 be an integer, 7 > 1, a > 
7 + - and 5 e (0,1). If u G H S ;](T x R+), then for any k = 1, 2, • • • , s, 



(A.5) 



[i + vV9k\\» < ll(i + vV&MW + r 2 ||(i + y y- l d k x {u - u)\\ 



L 2 



where gu := d x oo —d x (u — U). In addition, if u\ y=0 = 0, then for any k — 1, 2, • • • , s, 

to 

(A.6) ||(i + y yd k x co\W + ll(i +yy- 1 dt(u-u)\\<c 7 ^\\(i+ y yg4 L , + \\d k ML^} 

where C 1>a ^ is a constant depending on 7, a and S only. 

Assuming lemma [A2l which will be shown at the end of this appendix, for the moment, 
we can show inequality (A.4) as follows. 
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Applying lemma lA~2l for k = s, we obtain, from (A. 5) and (A. 6) 
1 

(A.7) 



2^ 4 ll(i + yrasWh <||(i + yVd°M\l 2 + ||(i + - £/)||i 2 



<C w {||(l + y)^||| 2 + ||^||i 2(T) }. 



Adding ||(1 + y y +a2 D a u\\ 2 L 2 to [(A~7)] we have 



[a|<s 

«1<S- 1 

^IkH^., < \\u\\ a E .„ + ||(1 + yy^u - U)\\ 2 L2 < C lja , 5 {\\uf H ^ + \\d s x U\\h { r)} 



which implies inequality (A. 4) 



Finally, in order to complete the justification of the almost equivalence relation (A.l) 
will prove the lemma IA.2I as follows. 



□ 



we 



Proof of lemma lATB . To prove (A. 5) , let us first recall from the definition of H*'g that 

6{l + y)- a <u<5-\l+v)-° 



(A. 



dyu<6-\l+y)- 



8yU 



so — — < 5 2 (1 + y) , and hence, for any k = 1, 2, 



UJ 



+ vV9h\\v< ll(i + yVd k M\ L , + ||(1 + y) 



d v oj , 



UJ 



+ y y- 1 d k x (u-u)\\ 



L 2 



< ll(i + vVd k M\ L . + r 2 || (i + y y- l d k x (u - u) \\ 



L 2 



which is inequality (A. 5) 

Next, 

9k ■= ud: 



Next, we are going to show (A. 6) The main observation is that we can also rewrite 

d k x ( u -uy 



UJ 



Ki+yy-^iu-u^KS-'iKi+yy 



. Thus, applying (A.8) t and part (ii) of lemma |RT| we have 

■U 2 

d k (u - U) 



uj 



d k U 



+ 



L 2 (T) 



(i + yy- a d y 



UJ 



L 2 



(A.9) 



< a h j- 2 {\\d k u\\L 2 m + ||(i + yyg k \\ L 2 }. 



Now, using triangle inequality, (A.8) and (A.9) , we also have 

ll(i + y yd k uj\\ L 2 < ii(i + yyg k \\ L 2 + r 2 n(i + yy- x d k x {u - u)\\» 
(a.io) < c w {||(i + y y gk \\ L2 + \\d k u\\ LHT) }. 

we prove (A. 6) 



Summing up (A.9) and (A.IO) 
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□ 

Appendix B. Calculus Inequalities 

In this appendix we will introduce several calculus inequalities for the incompressible 
velocity field (u,v), the vorticity u and the quantity g^. These inequalities are not related 
to any equations; they hold just because of elementary computations. 

B.l. Basic Inequalities. In this subsection we will state without proof two elementary 
inequalities (i.e., lemma [BTTI and lemma IBT21 below) . 
Let us first introduce Hardy's type inequalities. 

Lemma B.l (Hardy's Type Inequalities). Let f : T x R + — > R. Then 

(i) if A > — and lim f(x,y) = 0, then 

2 y— >+oo 

(B.l) ||(l +y )A / || i2(TxM+) <_^||( 1 + y) A+l^ / || i2(TxM+) . 



(ii) if A < — then 



(B.2) ||(l + y) A /ll^(TxR + ) < ^^ll/UI|L 2 (T)-^^||(l+y) A+1 ^/|U 2( TxR + ). 

The proof of lemma IB.ll is elementary, so we leave it for the reader. 
Next, we will state the following Sobolev's type inequality. 

Lemma B.2. Let f : T X 1R + — > KL Then there exists a universal constant C > such that 

(B.3) ||/||l°°(TxR+) < C{||/||l 2 (TxR+) + \\d x f\\L 2 (TxR+) + || ^/IU 2 (TxR+) } 

To prove lemma IB.2[ one may extend the domain of / to M 2 via the standard extension 



argument. Then inequalities (B.3) follows easily by the Fourier's inversion formula. We leave 
this to the reader as well. 

B.2. Estimates for H^] Functions. In this subsection we will use the weighted norm 
|| • ||h-|'T to control certain L 2 and L°° norms of u,v,u,gk and their derivatives. To derive 
these estimates, we shall apply lemma IB.ll and lemma IB.2| which was introduced previously 
in subsection IB.ll 

Our aim is to prove the following 

Proposition B.3 (L 2 and L°° Controls on u,v,cu and g^). Let the vector field (u,v) defined 
on T x M + satisfy the incompressibility condition d x u + d y v = 0, the Dirichlet boundary 
condition u\ y=Q = f L=o = and lim u = U . If the vorticity u := d y u G H^J for some 

constants s>4,7>l,a>7 + - and 5 G (0, 1), then we have the following estimates: there 
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exists a constant C S)7iCT ,s > such that 

Weighted L 2 Estimates: 

(i) for all k = 0, 1, • • • , s, 

(B.4) ||(1 + yy-'d^u - U)\\ L2 < C, lW {||u;||^ + II^IIW, 

(ii) for all k = 0, 1, ■ • • , s — 1, 
d k x v + yd k x +1 U 



1 + y 



(B.5) 

(iii) for all \a\ < s, 

(B.6) \\(l+ y y +a2 D a u\\ L 2 < 

(iv) for all k = 1, 2, • • • , s, 
(B.7) \\(l + yVgk\\v < 



<C s ^ s {H\ H ^ + \\d s x U\\ L ,}, 



i? 



(s,0). 



where the quantity := <9 x .u; —d x (u — U). 

to 



Weighted L°° Estimates: 

(v) /or a// = 0, 1, ■ • • , s — 1, 

(B.8) \\%u\\l<*> <C7 a , w {||w|| flS .T + 

(vi) /or a// = 0, 1, • • • , s — 2, 



CiwAMbF + II^IIW */ k = 1, 2, • • • , s - 1 

IMIhj.7 ifk = s, 



L 2 



}• 



<9£r 



<C S)7jCT>a {||a;|U r + ||^C/|U 2 }, 



(B.9) 

(vii) /or a// |a| < s — 2, 
(B.10) 

Proo/. 

(i) It follows from the definition of || ■ that ||(1 + y) 7_1 ^(w — U)\\ L 2 < \\u — [7|| HS , 7 -i. 



1 + y 

||(1 + y)-r+°*D a U>\\ Lao <C s Ju\\ H sn. 



so inequality (B.4) is a direct consequence of the almost equivalence inequality (A.l) 



(ii) Applying part (ii) of lemma iBTTl and inequality (B.4), we have 
d k v + yd k+1 U 



1 + y 



< 2\\d k x +1 (u - U)\\» < C^AM*" + \\%U\\v} 



L 2 
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which is inequality (B.5) 



(iii) Inequality (B.6) follows directly from the definition of ||w Hfl-s.-r and inequality (A.l) 

d TV 

(iv) Since u G H*'], we know that ||(1 + y)— — ||l°° < 8~ 2 , so using triangle inequality, 
inequalities (B.4) and (B.6) , we have 



which is inequality (B.7) for k — 1,2, ■ ■ ■ ,s — 1. When k — s, the better upper bound in 



(B.7) follows directly from the definition of ||o;||^.-> 



(v) For any k = 1,2, ■ ■ ■ , s — 1, applying lemma |B~2} inequalities (B.4) and (B.6) we have 



\d« x {u - U)\\ L - < C{\\d k (u - U)\\ L2 + \\d k+ \u - U)\\ L2 + \\d" x d y u\\ L >} 
<C s ^s{\\u\\ m „ + \\d s x U\\ L ,}, 



and hence, by the triangle inequality and ||c?*Z7||£,oo < C 8 ||^Z7||ia, we justify (B.8^ 



For the case k — 0, let us first recall from the hypothesis that u := d y u > 0, so 

" + 00 



(B.ll) 



< u < U 



u dy. 



Thus, using Cauchy-Schwarz's inequality and estimate (B.6) , we have 

1 



U 



L 2 



to dy 



dx < 



[i + v)M 



2 7 - 1 IIV ' *' " uL2 - 2 7 - 1 
and hence, by (B.ll), Sobolev inequality and H^C/H^ < (7 S ||9*?7||£^, we obtain 

\\u\\ L ™ < \\U\\ L °° < C{\\U\\ L 2 + \\d x U\\ L 2} < C 8 „ t<TtS {\\u Wh-,1 + \\d s x U\\ L 2} 



< 



\U}\ 



which is inequality (B.8) 



(vi) Applying triangle inequality, lemma IB~2| d x u + d y v = and u = d y u, we have 



1 + y 



< 



yd k x +1 u 



l + y 



d k x v + yd k x +1 U 



l + y 



<c{||^ +1 f/|| L2(T) + ||^ +2 f/|| L2(T) + 



d k x v + yd k x +1 U 



l + y 



L 2 



l + y 



L 2 



+ 



d k v + yd k+1 U 



+ 



L 2 



d k+1 {u-U) 



L 2 



d k+1 uo 



l + y 



L 2 



which implies inequality (B.9) because of Wirtinger's inequality and (B.4) - (B.6) 



(vii) Inequality (B.10) follows directly from lemma IBT21 and inequality (B.6) 



□ 

B.3. Estimates for Functions Vanishing at Infinity. In this subsection we will first 
define certain weighted norms involving time and spatial derivatives. Then we will state two 
basic inequalities about these norms, see proposition IB.5I below. Finally, in proposition IB. 61 
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we will control the weighted norms of u and v by that of u provided that u — U and its 
derivatives vanish at y — +00. The vanishing hypotheses (i.e., the decay rates) are usually 
guaranteed by proposition IC.ll in applications. 
Let us begin by defining the weighted norms. 

Definition B.4. For any s'gN and 7 6 1, we define 

It! [-1 

IN . |||2 ._ iioi . 112 IN . 1112 ._ II oi . 112 

III llls',7 ~ / j W U t Hifs'-2!,7(TxR+)' III \Ws' ■— /_^\\ u t llflV-2I(T)> 



1=0 



ll|-HI^7 : =E E lia+^^-Hi-M-H) and 

/=0 |a|<s-2Z 

HI . m|2 ._ || a . m2 

HI I lis', 00 -~ / ; \\ u t lliy s '- 2i .°°(T) 
1=0 

s' s' 
where [— ] denotes the largest integer which is less than or equal to — . 

Using Holder inequality and Sobolev inequality, one can easily show the following 

Proposition B.5 (Basic Inequalities for Weighted Norms). 

(i) For any s'eN and 7, 71, 72 € R with 7 = 71 + 72, 

II 1-^1-^2 1 ||s' ,7 < Cs'|||-^ 1 l|||s',oo,7i|||-^2|||s',72- 

(ii) For any s' > 5 and 7, 71, 72 G f with 7 = 71 + 72, 

|| I-F1-F2I ||s',7 < Cs',7i,7 2 {ll I -^1 1 l|s'-l,7i || 1-^2 1 ||s' ,72 + II I F\ 1 II s',71 1| | F"2. | || s'-l/ya } • 

Next, we will state the weighted controls on u and v as follows. 

Proposition B.6 (Weighted Controls on u and v). For any s' > 4 and 7 > 1, let the vector 

field (u,v) defined on T x M + satisfying the incompressibility condition d x u + d y v = 0, the 

s' 

Dirichlet boundary condition v\ y=0 = and lim ctd^u = d\d^U for all I = 0, 1, ■ • • , [— ] 

and k = 0, 1, • • • , s' — 21 + 1. Denote the vorticity u := d y u. Then there exists a universal 
constant C > such that 



I s',7 - 



(B.12) || \u — U\ \\ s 'fl < C|| |u;| || s / )7 and \\\v + yd x U\\\ s i-i < C\\\d x u\ 

Outline of the proof. The hypotheses of proposition IB .61 allow us to apply lemma IB. II to 
d\dl(u - U) and d l t d x (v + yd x U) provided that 21 + k < s', so we obtain |||w - U\\\ 8 r t0 < 
C|| \cj\ ll^i and || \v + yd x U\ \\ s '-i < C\\\d x co\\\ s ' t i which imply (B.12) since 7 > 1. □ 
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Appendix C. Decay Rates for H^'] Functions 

The aim of this appendix is to prove that the actual pointwise decay rates of H s a 'J functions 
at y = +00 are better than the decay rates obtained by the Sobolev embeddings. The 
proof relies on a pointwise interpolation argument (see lemma IC.3I below) , which is a direct 
consequence of the Taylor's series expansion. 

More specifically, we will prove the decay property of D a uj as y goes to +00 as follows. 



Proposition C.l (Decay Rates for H^] Functions). Let s' > 4 be an integer, 7 > l,a > 



7+2 an ^ ^ (0, 1). If ^ E H s a ^ in , then oj is s' + 2 times differentiable and there exists a 
constant CV, 7 ,5,||w|| s , +47 > such that for all \a\ < s' + 2, 

(C.l) \D a u\ < C s ,, iAMhs , + ^ (1 + vY ba inTx 

where the exponent 



(C.2) 



cr + (2H- 2 -i) 7 



2 |a|-2 



7 + a2 



if\a\ < 2 
+ a 2 if 2 < \a\ < s' + 1 
if \ a \ = s' + 2. 



Remark C.2 (Decay Rates from Sobolev Embeddings). Using the standard Sobolev embed- 
ding argument and the definition of H s aS ' 7 , one may prove that if u G H s ^\ then oj is 
s' + 2 times differentiable and there exists a constant C sn > such that 



(C.3) 



\D a u\ < 



6-\l + y)-°- 



C S)7 \\co\ 



H s '+ 4 <i 



-7— «2 



if |a| < 2 

if 2 < Id < s' + 2. 



Thus, the interesting part of proposition IC.ll is the decay rate of (C.l) is better than that 



of (C.3) This slightly better pointwise decay will help us to deal with the boundary terms 



at y = +00 while we are integrating by parts in the ^/-direction (cf remark I5T 



Proof of proposition I C. 1\ According to remark 10.21 we are only required to justify the 



inequality (C.l) with the decay rate defined in (C.2 



First of all, let us state without proof the following calculus lemma. 



Lemma C.3 (Pointwise Interpolation). Let f : TxIR + — > R be a twice differentiable function. 
Then we have the following: 

(i) If there exist constants Co, C 2, b andb 2 such that \d l x f\ < Ci(l + y)~ bi for alii = 0, 2, 
then 

b + b 2 



\d x f\<2^Co~C~ 2 (l + y y 



in T x 
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(ii) If there exist non-negative constants Co, C 2 , bo and b 2 such that \d y f\ < Cj(l + y) hi 
for all i = 0, 2, then 

\d y f\<2^C 2 (l + y)- bJ ^ mTxf+. 

The proof of lemma ICT31 is based on the standard Taylor's series expansion technique, and 
will be omitted here. 

Now, applying lemma ICT31 to D a u inductively on \a\ = 3, 4, ■ ■ • ,s' + l with the inequality 
(C.3), we prove (C.l) with the exponent b a defined in (C.2) □ 



Remark C.4 (Further Improvement on the Decay Rate). The decay rate b a defined in (C.2) 
is obviously not optimal because one can apply the pointwise interpolation lemma ICT31 again 
to further improve it. However, we do not intend to optimize it here. Indeed, repeatedly 
applying the pointwise interpolation lemma 10.31 one may improve the decay rate b a as 



cr + a 2 

(s' + 2- \a\)a+(\a\ - 2)7 



if \a\ < 2 
+ if 3 < |a| < s' + 1 
if \a\ = s' + 2. 



^7 + «2 

We leave this proof for the interested reader. 

Appendix D. Equations for a e and g\ 

In this appendix we will derive the evolution equations for a e 
d s x u e - a e d s x (u e - U) provided that to 6 > 0, (u e , v e , co e ) and (p £ , U) satisfy |(4.1)| - 1(4.4)[ These 
derivations just follow from direct computations. 



and g e s 



Equation for a e : 

Differentiating the vorticity equation (4.3) L with respect to y once, we obtain 
(D.l) (d t + u e d x + v e dy)d y u e = e 2 d 2 x d y u e + d*uj e - u e d x u e + d x u e d y u e . 



Using (D.l) and the vorticity equation (4.3) L , we can compute 

(d t + u e d x + v e d y )d y u e d y uj e (d t + u e d x + v e d y )uj € 



(d t + u e d x + v e d y )a € 



CO 



c\2 



(D.2) 

On the other hand, by direct differentiations only, one may check that 



<9V 



— d x u e + a e d x u e 



= d^yUf _ ^ _ ^ 



(D.3) 



Co> e 

dfa e = — — 



— a 



- 2a € d,a e . 
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(D.4) (ft + u e d x + v e d y - e 2 d 2 x - d 2 )a € = 2e 2 ^^d x a e + 2a e d y a £ -g[ + a € d x U, 



where g\ := d x u e — a e d x (u e — U). 



Equation for g\: (Derivation of equation (5.31)) 



Differentiating the evolution equations for ui e and u e — U (i.e., equations (5.29)) with 
respect to x s times, we have 

' (d t + u e d x + v e d y - e 2 d 2 x - d 2 )d s x u e + d s yd y co e 



(D.5) 



j=0 V/ j=i Vj/ 



(d t + u e d x + v £ d y - e 2 d 2 x - d 2 y )d s x {u e -U) + d s x v e u e 



= - E 0) % - to - E ( ■) ^TVav 

j=o \>J/ j=i Vj/ 

-E(-)^ ( "' '" r - 

3=0 W 

To eliminate the problematic term d x v € , we subtract a € x (D.5)k from (D.5) L , and obtain 
{dt + u £ d x + v e d y - e 2 8 2 x - d 2 y )gt + {(d t + u e d x + v e d y - e 2 d 2 x - d 2 )a e }d s x {u e - U) 

= 2e 2 d s x +1 (u e - U)d x a e + 2d s x u e d y a e - ("Xl'j ■ \<K V 
(D.6) j=0 \JJ 

-E (*Wv{^ e -a e ^ e } + a e £ ( S )di(u e -U)d x -i +1 U. 

1=1 7=0 



Substituting (D.4) into (D.6), we obtain equation (5.31) 



Appendix E. Classical Maximum Principles 
The main purpose of this appendix is to state two classical maximum principles, which 



are useful in subsection I5.2[ for parabolic equations. 

The first lemma is the maximum principle for bounded solutions to parabolic equations. 

Lemma E.l (Maximum Principle for Parabolic Equations). Let e > 0. If H G 

C([0,T];C 2 (T x R+) n C^QO, T]\ C°(T x E + )) is a bounded function which satisfies the 
differential inequality: 



{d t + hd X + b 2 dy - t 2 dl ~ d 2 }H < fH 



in [0, T] x T x 
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where the coefficients 61,62 and f are continuous and satisfy 
62 



1 + y 



(E.l) 

then for any t G [0, T], 
(E.2) 



< +00 



and 



L°°([0,T]xTxR+) 



L°°([0,T]xTxR+) < A, 



sup H(t) < max{e A '||#(0)|| L oo (TxM+) , max{e A( * r) \\H(T)\ y=0 \\ L ^ iT) }}. 

TxR+ re[0,t] 



The proof of lemma IE. 1 1 is a direct application of the classical maximum principle. For 
the reader's convenience, we will outline its proof as follows. 



Outline of the proof of lemma\KJ\ For any fx > 0, let us define "H := e xt H—Li 
/iln(l + y). Then one may check that for any t G (0, T], 



1 + y 



{d t + hd x + b 2 d y - e 2 d 2 - d 2 + (A -f)}H< 



in [0,t] x T x 



so by the classical maximum principle for parabolic equations, we have 

max U < max{ 11^(0)1^00^x15+), max{e~ AT ||if (t)|^ =0 ||l°°(t)}} 

[0,t]xTx[0,fi] re[0,t] 

provided that R > exp ( — \\H\\loo J — 1. Therefore, for any (x, y) G T x IR + , we have 

\» J 



H(t,x,y) - n 



1 + y 



e M t - Lie M ln(l + y) 



< max{e A *|| J ff(0)|| L ^ (TxM+) , max{e A( * r) \\H(T)\ y=0 \\ L ^ m }} 

rg[0,t] 



which implies (E.2) if we take the limit fi — > + and replace the arbitrary time t by t. 



□ 



The second lemma is the lower bound estimate on bounded solutions for parabolic 
equations. 

Lemma E.2 (Minimum Principle for Parabolic Equations). Let e > 0. If H G 

C([0,T];C 2 (T x R + )) n ^([O.rjjC^T x R + )) is a bounded function with 



and satisfies 



K(t) := mini min H(0), min H\ v=0 \ > 

TxR+ [0,t]xT y 



{9 4 + 61^ + b 2 d y - e 2 d 2 - d 2 }H = fH 



where the coefficients 61,62 and f are continuous and satisfy (E.l), then for any t G [0,T] ; 



(E.3) 



min Hit) > (1 - \te xt )K{t). 

TxR+ 
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The proof of lemma IE.2I is also standard and very similar to that of lemma IE. II We will 
outline it here for the reader's convenience as well. 



Outline of the proof of lemma 
e- xt {H-n(t)} + \\K(t)+ft 



KM For any fixed t G [0, T] and fi > 0, let us define h 
b 2 ) 

t + jj, ln(l + y). Then one may check that 



1 + y 



{d t + hd x + b 2 d y - e 2 d 2 x - dy + (\ - f)}h > in [0, ijxTx R+ 

so by the classical maximum principle for parabolic equations, we have 

min h > 

[0,t]xTx[0,i?J 

provided that R > exp ^— {||i?||ioo + — 1. Taking the limit R — > +oo, and then 

H — > + , we obtain 

H{t) > (1 - Xte xl )K{t) 



which implies inequality (E.3) if we replace the arbitrary time t by t. 

□ 
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